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We recall the concept of Baxterisation of an R-matrix, or of a monodromy matrix, 
which corresponds to build, from one point in the i?-matrix parameter space, the alge- 
braic variety where the spectral parameter(s) live. We show that the Baxterisation, which 
amounts to studying the iteration of a birational transformation, is a "win-win" strategy: 
it enables to discard efficiently the non-integrable situations, focusing directly on the two 
interesting cases where the algebraic varieties are of the so-called "general type" (finite 
order iteration) or are Abelian varieties (infinite order iteration). We emphasize the 
heuristic example of the sixteen vertex model and provide a complete description of the 
finite order iterations situations for the Baxter model. We show that the Baxterisation 
procedure can be introduced in much larger frameworks where the existence of some 
underlying Yang-Baxter structure is not used: we Baxterise L-operators, local quantum 
Lax matrices, and quantum Hamiltonians. 
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I. INTRODUCTION 

The Yang-Baxter equations are known to be a sufficient condition for the commutation of transfer 
matrices. Moreover, it has been shown that the commutation of transfer matrices necessarily yields 
a parameterization of the R-matrices in term of algebraic varieties ||, and that the set of inversion 
relations, combined together with the geometric symmetries of the lattice, yield a set of birational 
symmetries of the parameter space of the model, which are discrete symmetries of the Yang-Baxter 
equations HQ. Generically this set of birational symmetries is an infinite set. Combining these facts, 
one gets the following result: the Yang-Baxter integrability is necessarily parameterized in term of 
algebraic varieties having a set of discrete birational symmetries. An algebraic variety with an infinite 
set of (birational) symmetries cannot be an algebraic variety of the so-called "general type" |2j . In the 
following this set of discrete birational symmetries will be mainly seen as generated by the iteration of 
a birational transformation, thus canonically associating a discrete dynamical system. This birational 
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transformation can be of finite order, yielding a finite set of discrete symmetries : many Yang-Baxter 
integrable models correspond to this finite order situation || (see below : RSOS models ||, integrable 
chiral Potts model 0, free-fermion models, tetrahedron relations p|,p|,^0[ , ...). However, such birational 
transformations are generically (seen as a discrete dynamical system) infinite order transformations. 
The iteration of one point under such an infinite order birational transformation, yields an infinite 
number of points which can actually "density" an algebraic variety (elliptic curve, Abelian surface, 
... sec p"l] , |l2"| ] ) . One can thus deduce the algebraic variety from the iteration of one point by this 
birational transformation. This iteration procedure actually solves the so-called "Baxterisation" 
problem introduced by V. Jones in a framework of knot theory. The Baxterisation problem corresponds 
to actually find, from one (isolated) i?-matrix, satisfying the Yang-Baxter equation, a whole family of 
i?-matrices depending on one, or several, "spectral parameter(s)" satisfying the Yang-Baxter equation. 
In other words, the Baxterisation problem corresponds to actually build, from one point in the R- 
matrix parameter space (one Yang-Baxter integrable i?-matrix) , the algebraic variety where the spectral 
parameter(s) live. All this is also true for higher dimensional generalizations (tetrahedron relations 
...) of the Yang-Baxter equations, and for Baxterisation problems in dimensions greater than two |lc 

This paper will be illustrated by many examples of Baxterisation, in particular, the heuristic example 
of the Baxterisation of the sixteen vertex model which is, generically, non Yang-Baxter integrable, 
with a detailed analysis of the finite order situations of the Baxter model. One will finally show that 
the Baxterisation procedure can be introduced in much larger frameworks where the existence of some 
underlying Yang-Baxter structure is not clear, and not used. For instance, we will give several examples of 
Baxterisation of quantum Hamiltonians. We will also give several examples of Baxterisation of differential 
operators, starting with the simple example of the Baxterisation of the Toda L-operator, and then giving 
examples of Baxterisation of other simple local quantum Lax matrices. 



II. THE DISCRETE SYMMETRY GROUP ASSOCIATED WITH THE BAXTERISATION 

PROCEDURE 

Let us consider a quite general vertex model where one direction, denoted as direction (1), is singled 
out. Pictorially this can be represented as follows: 



(1) 

(1) 

where i and k (corresponding to direction (1)) can take q values, while J and L, in the other direction, 
take m values. 

One can define a "partial" transposition on direction (I) denoted t\. The action of t\ on the i?-matrix 
is given by @||: 

(hRf k i = R\l (2) 
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The i?-matrix is a (q m) x (q m) matrix which can be seen as q 2 blocks which are m x m matrices: 



R = 


/4M] 

42,1] 
A[3,l] 


41,2] 
42,2] 
43,2] 


41,3] • 
42,3] • 
43,3] ■ 


•• A[l,q]\ 

•• 42, g] 
■■ 43, g] 




\A[q,l] 


A[q,2] 


49, 3] • 


•• 4g,g]/ 



(3) 



where A[I,1], A[I,2], A[q,q] are m x m matrices. With these notations the partial transposition 
ti amounts to permuting matrices 4a, /3] an d 4A a ]> f° r au t nese block matrices. We use the same 



2 



notations as in Jl 4|J 1 5| Jl 6[ | , that is, we introduce the two following transformations on matrix R, the 
matrix inverse I and the homogeneous matrix inverse /: 

T:R — ► R- 1 , I:R — > det(i?) • iT 1 (4) 

and we introduce the (generically) infinite order homogeneous, and inhomogeneous, transformations : 

K = h-i, k = ti-T, (5) 

The homogeneous inverse / is a polynomial transformation on each of the entries my of R, which 
associates to each my its corresponding cofactor. Transformation / is an involution, whereas I 2 = 
(det(R)) qm _ 2 -X, where X denotes the identity transformation. Transformation t\ is also an involution. 
The transformation K is clearly a birational transformation on the entries my since its inverse transfor- 
mation is I ■ t which is obviously also a rational transformation. Transformation K is a homogeneous 
polynomial transformation on the entries my of the i?-matrix. 

For such vertex models of lattice statistical mechanics, transformations I and t\ come from the 
inversion relations jl7 18 1, and the geometrical symmetries of the lattice, in the framework of integrability 



and beyond integrability. These involutions generate a discrete group of (birational) automorphisms of 
the Yang-Baxter equations ||,|| and their higher dimensional generalizations [HJ. They also generate 
a discrete group of (birational) automorphisms of the algebraic varieties canonically associated with 
the Yang-Baxter equations (or their higher dimensional generalizations) Q]. In the generic case where 
the birational transformation K = t\ ■ I is an infinite order transformation, this set of birational 
automorphisms corresponds, essentially, to the iteration of K. When the infinite order transformation 
K densities the algebraic varieties |l^] , one can deduce the equations of these algebraic varieties exactly 
and thus solve the so-called " Baxterisation problem" . 

Terminology : In the following, we will call "i?-matrix" the matrix we iterate with (^), even if 
this "i?-matrix" does not satisfy a Yang-Baxter relation. From now, on we will call "Baxterisation" 
the procedure which amounts, for a given i?-matrix, to iterating transformation (^), and finding the 
algebraic variety which contains all the points of this iteration. A good example corresponds, for instance, 
to get, from isolated i?-matrices satisfying the tetrahedron relations |T3], i?-matrices also solutions of 
the tetrahedron relations but now depending on spectral parameters. 



III. BAXTERISATION : A WIN- WIN STRATEGY 



When one "Baxterises" an R-matrix one can get three different kinds of situations : either the orbits of 
(||) are chaotic, and therefore one cannot expect any "nice" parametrization of the lattice model, or the 
orbits correspond to algebraic varieties, and one can actually introduce some "well-suited" parametriza- 
tion of the model associated to these algebraic varieties, which will be extremely precious for any further 
analysis of the model (analysis of the Yang-Baxter equations, or higher dimensional generalizations, 
calculations of partition function per site, ...), or, finally, these orbits are finite orbits. When the itera- 
tion of (^) is infinite and yields algebraic varieties, one can show that the algebraic varieties are not of 
the so-called "general-type" [|| (using the algebraic geometry terminology). For instance, for algebraic 
surfaces, one can actually classify the various surfaces that are not of the "general type" : product of 
elliptic curves, Enriques surfaces, Kummer surfaces, Abelian surfaces, ... When the orbits are finite one 
is back to algebraic varieties of the "general-type" , which is a very large set of different situations, quite 
hard to classify However, writing for some integer N, the (projective) condition : 

K N (R) = C-R (6) 

actually gives the equations of the algebraic varieties corresponding to this very condition (^). For 
instance, one can find the algebraic subvariety of the chiral Potts model on which the higher genus 
Au- Yang-Baxter-Perk solutions p0| live, as a finite order condition M. 
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To sum up : one has the following situations : 

• Either the group is infinite, and one gets the "precious" parametrization of the model (even if it is 
a parametrization in terms of theta functions of g variables, ...). 

• The group is finite, and one actually gets the equations of the algebraic varieties corresponding to 
this situation, writing condition (^). This finite order situation is a "gold mine" as far as integrable 
situations are concerned HQ. Let us recall, for instance, the example of the tetrahedron solutions H|). 

• The group is infinity, but the orbits, corresponding to the iteration of (^), arc chaotic : one does not 
have algebraic varieties El and therefore one does not have any Yang-Baxter structure 0. This situation 
corresponds to an exponential growth of the iteration calculations [^1 22 , 2^] . 



In every case the Baxtcrisation procedure helps to avoid the points of the parameter space, corre- 
sponding to exponential growth of the calculations, where no Yang-Baxter relation can be expected, and 
enable s to get analytically the algebraic varieties corresponding to finite order iterations (see section 



( VIA ) below), or enables to actually build the algebraic varieties from the iteration of an infinite order 



birational transformation. 

The Baxterisation procedure amounts to studying the iteration of a (resp. several) birational (resp. 
polynomial) mapping which actually corresponds to discrete symmetries of the parameter space of the 
models. It thus provides a natural link between lattice statistical mechanics (field theory, ...) and the 
theory of discrete dynamical systems. Furthermore, it also provides links with many other domains 
of mathematical physics, or mathematics. For instance, as far as effective algebraic geometry, or even 
arithmetic, is concerned, the Baxterisation provides many natural, and simple, examples of Abelian 
varieties with an infinite set of rational points p"^ . 

Let us first recall some well-known Yang-Baxter integrable situations corresponding to finite order 
conditions. 



A. Finite order conditions : free-fermion conditions 

The matrix R of the asymmetric eight- vertex model is of the form : 

/a d'\ 

\d a'J 

The free-fermion condition is : 

aa' - dd' + W -cc' = (8) 
A matrix of the form (Q) may be brought, by similarity transformations, to a block-diagonal form: 

R= ( R i? 2 )> With Rl= (d a') and R *=( b c b> 

If one denotes by 8\ = aa' — dd' , and by 82 = W — cc' , the determinants of the two blocks, then the 
homogeneous inverse I (polynomial transformation) just reads : 



2 One may have algebraic subvarieties compatible with this situation, the subvarieties corresponding to the two 
previous situations. 
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One can, however, try to characterize the "complexity" of t 



entropy, or the Arnold complexity, of the birational mapping (|5) one iterates [ pljp^Jj3| 



lis chaotic situation, calculating the topological 
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a — ► a' ■ §2, a' — ► a ■ S2, d — > — d • (52, d' — * — d' • #2, (9) 
6^6' hf->b-6i, c^-c-Si, d -»■ -c' • Si 

and transformation ii is given by : t\ : c <-> d! and : d <-> cf. The condition (^) is left invariant by ti, 
/, and thus K = t\ ■ I. It is straightforward to see that condition (||) is Si = —82 and has the effect 
of linearizing I into : 

a^a', a' -ki, d-*--d, d'-^-d', b->-b', b' -> -b, c^c, c' -> c' (10) 

The group, generated by / and t\, is then realized by permutations of the entries, mixed with changes 
of signs, and its orbits are thus finite. The situation depicted here, namely an inversion relation that 
reduces, on some algebraic subvariety, to permutation of the entries up to signs, also occurs for free- 
fermion two-dimensional vertex models on a triangular lattice (see Sacco and Wu in 0]), or for the 
three-dimensional vertex corresponding to the Zamolodchikov-Baxtcr solution of the tetrahedron equa- 
tions |p5| , p6| , p7| . Details can be found in j^]. 

The finite order situations are extremely favorable for Yang-Baxter integrability (chiral Potts model 0] , 
RSOS models ||, tetrahedron solution ||, free-fermion solutions, free para-fermions, ...). If one tries to 
find new Yang-Baxter integrable models, one should certainly first try to write, systematically, all the 
algebraic subvarieties corresponding to these (projective) finite order conditions. However, in this finite 
order case, one may say that the "Baxterisation procedure" does not work, or "works to well" : it is too 
degenerate. As a consequence, many remarkable results, and structures, occur (polynomial representa- 
tion of the natural i nteger s together with their multiplication, ...). These results, and structures, will be 



sketched in section ( VI A ) 



B. The group is infinite : iterations associated with the sixteen- vertex model 

Let us now consider, with the example of the sixteen vertex model Eft , the situation where K, or 
K, is infinite order, thus yielding a non-trivial Baxterisation. 

In the case of 4 x 4 matrices (see Fig but with q — m = 2), a particular permutation of the entries 
of the matrix, t\, has been introduced in the framework of the symmetries of the sixteen-vertex model 
p8[ . The action of two partial transpositions ti , and £2 , on the i?-matrix is given by [p8| : 

(*i*)S = 4 J \ faR)%=I$j, t = ti-t 2 (11) 

Transformation t is nothing but the matrix transposition : t commutes with the matrix inversion /. If 
one denotes the entries of the i?-matrix, this permutation corresponds to: 

ti : mi 3 <-* m 3 i, mu <-> m 32 , m 2 3 «-> ma, m 2 A <-> m 42 (12) 

which amounts to permuting the two 2x2 (off-diagonal) sub-matrices of the 4x4 i?-matrix. This 
transposition ti corresponds to a partial transposition of one direction (say the horizontal one denoted by 
"1" , the other transposition corresponding to the other direction denoted by "2" ) of a two-dimensional 
vertex model [|]J4|,|28| (see Fig.|l|). 

Remarkably, the symmetry group, generated by the matrix inverse / and transformation t\, or by 
the infinite generator K = ti • I, has been shown to yield elliptic curves [Q,^8| which foliate the whole 
parameter space of the sixteen vertex model. One should not confuse the integrability of the symmetries 
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of the parameter space of the sixteen vertex model (namely the mappings considered here) and. the 
Yang-Baxter integrability |28[| : the sixteen vertex model is not generically Yang-Baxter integrable H. 

The integrability of the birational mapping K, or, equivalently, of the homogeneous (bi-)polynomial 
transformation K, is closely related to the occurrence of remarkable factorization schemes |O,0. In 
order to see this, let us consider a 4 x 4 matrix Mo — R, and the successive matrices obtained by 
iteration of transformation K = t\ • J, where t\ is defined by (jl^)- Similarly to the factorizations 
described in one has, for arbitrary n, the following factorizations for the iterations of K : 



M _ K(M 1l+1 ) _ det(M n+1 ) ~ _ K(M n+2 ) _ M n+3 

M n+2 - 75 , Jn+2 - 75 , A iil^n+2 — , 7777 7 — ~? 7 (^J 

In In det(M„ +2 ) /n+l/n+3 

where the /„'s are homogeneous polynomials in the entries of the initial matrix Mq and the M ra 's are 
"reduced matrices" with homogeneous polynomial entries pl|^2[ . 

Let us denote by a n the degree of the determinant of matrix M n , and by (3 n the degree of polynomial 
/„, and let us introduce a{x), (3{x) which are the generating functions of these a„'s, /3 ra 's : 

oo oo 

a(x) = ]T a n ■ x n , (3(x) = ^/3 n ■ x n , (14) 

71 = 71 = 

From these factorizations, one can easily get linear relations on the exponents a n , f3 n and exact expres- 
sions for their generating functions and for the q„'s and f3 n 's: 

4fl + 3a; 2 ) 4x 

a ( x ) = "Ti vT> l 3 ( x )= 7i vT> ««= 4(2n 2 + l), f3 n = 2n(n+l) (15) 

(1 — xy (1 — xy* 

One has a whole hierarchy of recursions integrable, or compatible with integrability i jOj . For instance, 
one has : 

fn In+3 ~ /"+4 In+1 _ In+1 I n +4 ~ /w+5 /n+2 ^g-j 



/n-1 /n+3 /n+4 — /« In+1 /n+5 In In+4 /n+5 — In+1 /n+2 /n+6 

In an equivalent way, introducing the variable x n — det(K n (R)) ■ det(K n+1 (R)), one gets a hierarchy 
of recursions on the a;„ 's, (see |ll[), the simplest recursion reading : 

X n +2 — 1 X n +1 — 1 2 



^n+l a?7%+2 ^71+3 — 1 £n+l x 7i+2 — 1 



%n+2 (^) 



Equation ( |l7|) is equivalent to ( |l6| ) since x n = (/^ • fn+2)/(fn+i ' /n-i)- It can be seen that these 
recursions ( |16[) and (|l7| ) are integrable ones |pJ| , For this, one can introduce |l4[ a new (homogeneous) 
variable : 

In+1 ' /n— 1 ,i 9n+l /io\ 
?7i = J2 tn(2n 1 ^ = ( 18 -' 

In 9™ 
and end up, after some simplifications, with the following biquadratic relation between q n and q n +i ■ 

il-il+i + v-qn-qn+i + p- (q n + q n +i) - A = o (19) 



5 Using this elliptic parameterization one could imagine that the "inversion trick" together with some 

"well-suited" analytical assumptions, could allow to actually get the exact expression oftne partition function per 
site of the sixteen vertex model beyond the Yang-Baxter integrability. The model would be "calculable" without 
being Yang-Baxter integrable : this remains an open question. 
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which is clearly an integrable recursion [ fLl| . In terms of the /„'s, the three parameters p, A and /i 
read : 



fl 2 fi — /2 2 fz . /2/4— /l /3 2 f2 5 — f3 Z fl' i — fl'fi + fs fi" 



/' 



/f/3 - fi ' tf/s-ZI ' n hh-(fih-fi) 



It may also be interesting to introduce : 

4 • A + ^ 2 



(20) 



IV. B AXTERISATION OF MONODROMY MATRICES : 2 M x 2 M MATRICES 



Let us now recall the more general vertex model (see Fig. [j] above), where the singled out direction 
(1) corresponds to a two- dimensional "auxiliary space" (that is q = 2). The action of t\, the "partial" 
transposition on direction (1), is given by p8j : 

yij u kJ . 1 A B\ {AC 



(hR)\ J L = ^ that is : t! : I ^ ~ I — I £ ~ I (21) 



where A, B, C and D are m x m matrices. It is a straight calculation to see that the matrix inversion 
reads : 



A B\ _ ((A-B-D-i-Cy 1 (C-D-B-^A)' 1 



C D 



B — A - C ■ D [D — C ■ A ■ B 



(22) 



This general framework enables to take into account the analysis of N-site monodromy matrices 
(take m = 2 N ) of two-dimensional models, as well as the analysis of d-dimensional 2 d -state vertex models 
(take m = 2 d ~ 1 ). Let us just give here a pictorial representation of the two sites (N = 2) monodromy 
matrix of a two-dimensional model and of a three-dimensional vertex model: 




Denoting s — 2m the size of the matrices, the analysis of the corresponding factorizations yields for 

arbitrary n, "string-like" factorizations Jl^Jl^]. For arbitrary m (equal to 2 d_1 or not), the analysis of 
the factorizations of the iterations of transformation K yields: 

71 r K(M) f , f(M) f dct(A/ 1 ) KjAh) dot(M 2 ) K(M 2 ) 

Mi = K(M ), /1 = det(Alo), h = JS _ 4 , M 2 = JS _ 5 , h = -pf — 7^34, = 



fS—i ' z fs— 5 ' JO f7 rs~ 4' f§ rs 

Jl Jl Jl ' J2 Jl ' J2 



_ det(M 3 ) _ ir(M 3 ) det(M 4 



l»—*> . f7 . fs-4, ' j2(s-5) j 5 . j-s-5 ' J j?8 . f 2 . jr7 . fa-i 



and, for arbitrary n, the following "string-like" factorizations : 
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K(M n ) — M n+ i ■ f* 5 • f^_ x ■ f^ 2 5 ' • f®_ 3 • /^_f 4 5 • /^_ s • • • , 

det(M n ) = f n+ i ■ f^ 4 ■ f n _ x ■ f}^ 2 4) • f®_ 3 ■ /^_f 4 4) • /®_ 5 ■ 4) • • • (24) 



One easily gets from (BJ) : 



• s 2 2; (1 + a; 2 ) a , si 



= i — i — + s T\ — i — m \i' 



l + s (1 + x) (l-x) 4 ' ^ w (1-z) 3 



Q'„ = (2n+l)(2n 2 + 2n + 3) , Ai=|n(n + 1) (25) 

The a„'s and /3 ra 's are, respectively, cubic and quadratic functions of n. 



A. Towards Bethe Ansatz : the propagation property 



This polynomial growth of the calculation can be understood as follows. One of the "keys" to the 
Bethe Ansatz is the existence (see equations (B.10), (B.lla) in ^9|) of vectors which are pure tensox. 
products (of the form v ® w) and which R maps onto pure tensor product v' (£> w' . This key property □ 
was called propagation property b« R. J. Baxter, and corresponds to the existence of a Zamolodchikov 
bra ]3S| for the Baxter model 0. This "propagation" equation reads here : 

R ( u ® V) = n-u'®V' with : u= ( 1 ) , u = ( \) , (26) 



p I \p 



vectors V and V having m coordinates. One can rewrite ( |26|) under the form 

c d) { P V v) = vy/.v 1 (27) 

Actually, for all the vertex models for which transposition t\ can be represented as (^l]) (namely 
monodromy matrices, or d-dimensional vertex models, with "arrows" taking two colors, ... ), one can 
associate an algebraic curve of equation : 

det(Ap' - C - Dp+ pp' B) = (28) 

which form is invariant by t\ , / and thus by K or K 2 . As a byproduct this provides a canonical 
Jacobian variety for such vertex models, namely the Jacobian variety associated with curve (pq). This 
procedure, which associates with an i?-matrix the algebraic curve (|28|), originates from a key "propaga- 
tion" relation (p6|), closely related to the action of the birational transformations K p8|]. 



6 Which is "almost" a sufficient condition for the Yang-Baxter equation. In the case of the Baxter model this 
non trivial relation corresponds to some intertwining relation o f the p roduct of two theta functions, which is 
nothing but the quadratic Frobenius relations on theta functions |30lplJ. 

7 The existence of a Zamolodchikov algebra is, at first sight, a sufficient condition for the Yang-Baxter equations 
to be verified. Theta functions of g variables do satisfy quadratic Frobenius relations [pl[ , consequently yielding a 
Zamolodchikov algebra parameterized in term of theta functions of several variables. However this Zamolodchikov 
algebra is apparently not sufficient for Yang-Baxter equations to be satisfied |B3| . For the corresponding vertex 
models ( g replicas of the Baxter model coupled together) it may be possible thatThe associated i?-matrices, which 
are parametrised in term of theta functions of several variables (and thus correspond to a "nice" Baxterisation), do 
not satisfy the Yang-Baxter equations, but could be such that the partition function per site could be calculated 
exactly using the "inversion trick" : this remains an open question. 
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B. Continuous symmetries generalizing the gauge symmetries 



The birational transformation K = t\ - 1 can be represented as follows 



A B\ A' B 




-1 / v-l' 

B- A-C- 1 -D] 



= i c »)- (a „)- i>„ „ ; r \ D _ . A -,.^ (2<J) 

Let us introduce the following SL(m) x SL(m) transformation : 

aim), (a b\ (g^-a-g^ g^-b-g^) 

' \C D) [g^-C-G^ G^.D-G™ J 

where G^ and G^™' are two SL(m) matrices. It is straightforward to see that : 



G (m) . A , G (m) G (rn) . B ■ G^ \ / (G^)" 1 ' A ' ' (G^) _1 (G^)" 1 ' ' (G^) -1 

G M _ c G (m) G M _ D _ G (m) II (G M)-1 . c , . ( G (™))-i (G^)- 1 • D> ■ (G^)- 1 



and, in a second step, that : 



g {m \k 2 {R)) = k 2 {g {m) {R)) (31) 



In fact, such a result is not specific of a two-dimensional auxiliary space. Recalling the (qm) x (qm) 
i?-matrix (|3j), and the associated partial transposition ti, and introducing a SL(m) x SL(m) transfor- 
mation which transforms each to x to block A[a, f3] into Gl ■ A[a, 0] ■ Gr, one recovers, again, 
relation (flj). Of course there is nothing specific with t\, and one finds the same results for the partial 
transposition t% corresponding to the vertical line, with this time, a commutation between K 2 and a 
SL(q] x SL(q) transformation Q^ q > . Since £2 = ^i ■ t, and since t commutes with /, one immediately 
gets B with obvious notations the following SL(q) x SL(m) x SL(q) x SL(m) symmetry for K 2 : 



K 2 



(g [ Z ] ® G^ m) • i? • G^ } ® G^ m) ) = G [ l ] ® G^ m) • (if 2 (i?)) • G^ } <g> G^ m) (32) 



This symmetry drastically generalizes the well-known gauge symmetries p4j which correspond to 
G^ = (G^) -1 and G^ = (G^)^ 1 . Let us note that this kind of symmetries, generalizing 
the gauge transformations, have actually been used by R. J. Baxter to map the Baxter model onto an 
inhomogeneous six vertex model, in order to build the Bethe Ansatz of the Baxter model ]29|. 



V. SIXTEEN VERTEX MODEL AND BAXTER MODEL : REVISITING THE ELLIPTIC 

CURVES 



Considering the (non-generically Yang-Baxter integrable) sixteen vertex model, one finds that a canon- 
ical parameterization in terms of elliptic curves occurs in the sixteen homogeneous parameter space of 
the model pa] . This canonical parameterization is obtained □ from the Baxterisation procedure p8| . In 



8 Similarly for a three-dimensional 2 3 x2 3 i?-matrix one gets easily K 2 (g 1 ^(&g 2 ^®g 3 ^-R,3D-giR®g2R®g'iR) = 

9il ® 92L ® 93L ■ K 2 (R 3D ) ■ giR ® g2R ® g-sR. 

9 In the Yang-Baxter integrable subcase, the Baxter model, this elliptic parameterization, deduced from the 
Baxter isat ion procedure, is actually the elliptic parameterization introduced by R. J. Baxter to solve the Baxter 
model |B5|. 
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fact several elliptic curves ( assoc iated with different "spaces") occur : one corresponding to the factor- 
ization analysis of section ( III B| ), another from the propagation property (p6|), and another one from 
the iteration of the birational transformation K 2 in the sixteen homogeneous parameter space of the 
model p8| . Let us analyse the relations between these various elliptic curves and show that they actually 
identify. 

Let us recall the results and notations concerning the sixteen vertex model p8[ . Let us use the 
following notation for R : 



R 



I a\ a 2 bi b 2 

a 3 a 4 b 3 b 4 

ci c 2 di d 2 

\c 3 c 4 d 3 d it 



(33) 



A. Propagation property for the sixteen vertex model and the Baxter model 



Considering the sixteen vertex model (|33|), the propagation relation (|26|), for m = 2, becomes 
R (v n ® w n ) = /i • v n+ i ® w n+ i where : 



V n +1 



1 

Pn+1 



W n +l 



1 

Pn+1 



(34) 



and yields, by eliminations of p n , p n +i (resp. p n , p n +i) 1 the two biquadratic relations 

U + kl ■ Pn ~ h2 ■ Pn+1 + h ■ P n + k ■ P n +1 - ( l 9 + h») ' Pn ' Pn+1 

~ Il3 ■ Pn ■ Pn+1 + ho ■ Pn ' P n+1 + h " Pn ' Pn+1 = 



(35) 



h + he ■ Pn ~ h5 ■ Pn+1 + h ■ P n + h ■ Pn+1 ~ (^9 " hs) ' Pn ' P 



■ Pn+1 



ll7 ■ Pn ■ Pn+1 + hi ■ PnPn+l + l 6 P„ Pn+1 



K ^2 



^2 -2 







(36) 



where the U 's are quadratic expressions of the entries ( p3| ) of the i?-matrix p8[ . These two biquadratics 
can be seen as : 



P 
P 
1 



= 



and 



r g 2 



= o 



where the two 3x3 matrices read : 



4 X) = 



h ho h 

-I12 —{h + ^is) ~ '13 
h hi h 



and 



p(2) 

K 3 - 



I5 hi h 

~hh —{I9 — Us) — h7 
h he h 



(37) 



In the Baxter limit these two matrices reduce to a only one 3x3 matrix 



R 



bax 



J x H~ Jy J x Jy 

4J Z 

Jx Jy Jx ~t~ Jy 



(38) 



See the propagation property for the Baxter model, namely (B.10), (B.lla) in 
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where J x , J y , and J z are the three well-known quadratic expressions of the XYZ Hamiltonian : 
J x = a ■ b + c ■ a, Jy — a ■ b — c ■ a , J z = , 



Some i4T 2 -invariants can be deduced from SX(3) invariants of the two 3x3 matrices (37), namely a 

quadratic expression in the l^s (for for instance) h ■ h + h ' h — ho ■ — 1 12 ' '13 + ('9 + ^is) 2 > 
a cubic (in the U's ) which is nothing but their determinant, and a quartic one. Eighteen (algebraically 
related) quadratic polynomials (pi, . . . ,pis) which are linear combinaisons of the Zj's, and transform very 
simply under ti and /, have been found Introducing the ratio of these covariants pi's, one gets 

invariants of K 2 thus giving the equations of the elliptic curves : the elliptic curves are given by the 
intersection of fourteen quadrics Eg] . 



B. Reduction of a sixteen vertex model to a A' -effective Baxter model 

From (p3) for q = m = 2, one deduces a SI2 x s?2 x s^ 2 x s?2 symmetry on the .sixteen vertex 
model |2 ^^Furthermore, the i?-matrix of the sixteen vertex can actually be decomposed Eil as : 

Rsixteen = 9\h ® 92L ' ^Baxter ' 9lR ® 9lR (39) 

where RBaxter denotes the i?-matrix of an "effective" Baxter model and gm, g2R, 9il, 92L are 2x2 
matrices. The sixteen homogeneous parameters of the sixteen vertex are thus decomposed into four 
homogeneous parameters of an "effective" Baxter model and four times three parameters (four homo- 
geneous parameters) of the various 2x2 matrices : giR, g2R, giL, 32L- Using this very decomposition 
(p9|), and the previous symmetry relation ( |32| ) for q = m = 2, one actually gets : 

K 2 (R slxteen ) = gll ® g~l ■ K 2 (RBaxter) ■ giR ® g2R (40) 



The matrices g\R, g2R, giL, 52L of the decomposition ( p9| ) can thus be seen as constants of motion 
of the iteration of K 2 . 

If RBaxter belongs to a "special" manifold, or algebraic variety, Rsixteen, given by (|39|), will also belong 
to a "special" manifold, or algebraic variety : for instance, if Rsaxter belongs to a finite order algebraic 
variety for the iteration of K 2 , namely K 2N (RBaxter) = V ■ RBaxter, then Rsixteen will also belong to a 
finite order algebraic variety for the iteration of K 2 : K 2 N (Rsixteen) — V ' Rsixteen- If RBaxter belongs 
to a critical variety then Rsixteen given by (|39| ) should also belong to a critical variety. This last result 
does not come from the fact that g\h, g2L, giR, g2R are symmetries of the partition function (they are 
not, except in the gauge case : g±L — giR with 52L = 52i?) : they are symmetries of K 2 which is a 
symmetry of the critical manifolds. Therefore they are symmetries of the critical manifolds even if they 
are not symmetries of the partition function. 

A decomposition, like ([39|) , is closely associated to the parametrization of the sixteen vertex model in 
elliptic curves Q : given RBaxter, 9il, 92L, 9ir and g 2 R, one can easily deduce Rsixteen- Conversely, 
given Rsixteen, it is extremely difficult to get RBaxter, 9xl, 92L, ffifl^and g 2 R, however, and remarkably, 
it is quite simple to get RBaxter- Since giL, giR, g2L, g2R are /^-invariants, one can try to relate, 
directly, the "-ftT 2 -effective" covariants J x , J y and J z with the i4T 2 -invariants related to the recursion 
on the x n 's or the q n 's, namely p, fx, A, or n (see (pj|), (53)). In terms of these well-suited algebraic 
covariants, the previous parameters read : 



Finding, for a given R- matrix of the sixteen vertex model, the elements of this decomposition, namely Rb u 



and gi R , g 2 R, giL, g2L, is an extreme 



difficult process that will not be detailed here. Conversely, one can show 



easily that the matrices of the form (391) span the whole space of 4 x 4 matrices. 
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p= 4 J 2 J 2 J 2 , n= -2- {Jpl + J 2 z J 2 y + Jpl) , « = 4 • (J 2 Z + Jj + J 2 ) 

A = - ( ^ + J|j y 2 + J^Jy) 2 + 4 • + J 2 + ■ Jlj 2 x j 2 y (4l) 

One immediately recognizes some symmetric polynomials of J 2 , J 2 , J 2 . Therefore it is easy to see 
that the J Xl J y , J z can be straightforwardly obtained from a cubic polynomial P(u) : 

P(u) = 4:-u 3 -K-u 2 -2-(i-u-p = 4 • (u- J 2 ) ■ (u- J 2 ) ■ (u- J 2 ) (42) 

This is remarkable, because trying to get J x , J y , J z , by brute- force eliminations from (|39|), yields huge 
calculations. In fact, one only gets, from ([42]), the squares of the J x , J y , J z , but the critical manifold, 
as well as the finite order conditions (see below), only depend on J 2 , J 2 , J 2 . 



C. Biquadratic (19) versus biquadratic (35) 

The Baxterisation process is associated with the iteration of K, or, rather, K 2 . Since, as far as 
K 2 is concerned, one can reduce a sixteen vertex model to an "effective" Baxter model, one can try 
to revisit, directly, the relation between the biquadratic ( |l9| ) and the "propagation curve" ( |35| ) for the 
Baxter model. For the Baxter model relation ( |l9| ) becomes the biquadratic : 

3« ?n+l - 2 ■ (JlJl + JzJy + J 2 Jl) ' 9n ' 9n+i + 4J, 2 Jj J 2 • (q n + q n +i) 

+ (J 2 Jl + J 2 J 2 + J 2 J 2 ) 2 - 4 • (J 2 + J 2 + J 2 ) ■ J 2 J 2 J 2 = (43) 



which should be compared with the "propagation" biquadratics ( p5| ) of the Baxter model [|29|,|35| : 

TliPn, Pn+l) = {Jx - Jy) ' (pIpI+1 + 1) - (Jx + Jy) ' [Pn + Pn+l) + 4 ' J » ' Pn ' Pn+l 
= (Pi ~ 1) ' (pI+i -!)•■/«- (Pi + 1) ' (Pn+l + 1) • Jy + 4 • J z ■ Pn • Pn+l = (44) 

It is known that simple "propagation" curves, like (^i[), have the following elliptic parameteriza- 
tion |2||||: 

p n = sn(u n ,k), p n+ i = sn(u n+ i,k) where: u n+ i = u n ± A (45) 

where sn(u,k) denotes the elliptic sinus of modulus k and A denotes some "shift". The modulus |39] 
k is equivalent to the following modulus which has a very simple expression in terms of J 2 , J 2 and jf : 

J x J z 

At first sight it seems that one has two different elliptic curves (biquadratics), namely ( [43| ) which is 
symmetric under permutations of J 2 , J? and J 2 , and (^) which breaks this symmetry. Let us also 
consider the same "propagation curve" (|44|), but now between p n +i and p n +2, and let us eliminate p n +i 
between these two algebraic curves. One gets, after the factorization of (p n — p n +\) 2 ■ 

T 2 (pn, Pn+2) = (47) 
= 2 J 2 - (J 2 -J 2 )- (p 2 n p 2 n+2 + 1) + 2 J 2 J 2 ■ (p 2 n + p 2 l+2 ) + 4 ■ (J 2 J 2 - J 2 J 2 - J 2 J 2 ) ■ pn ■ Pn+2 

= (pi - 1) • (p 2 l+2 - 1) ■ 4 2 > - [pi + 1) ■ (p 2 l+2 + 1) ■ + 4 • J^ ■ Pn ■ Pn+2 - 

where J x 2 ^ , Jy 2 * 1 and J z 2 ^ are given below. One remarks that fl47| ) is actually of the same form as 
(p4]). The two biquadratic curves ( |43| ) and (^) are (birationally) equivalent, their shift A and modular 
invariant [p6[ being equal. Actually one can find directly the homographic transformation 
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a-Pn + (3 OL-Pn+2 + P , , 

In = : j , q n +i = — -j 48 

7 • Pn + 7 • Pn+2 + 

which maps ( f43| ) onto (|47|), the parameters a, ... S of the homographic transformation (48) being quite 
involved. The calculations are quite tedious and will be given elsewhere. 

The previously described elimination of p n +i, changing Ti into I^, amounts to eliminating u n+ i 
between u n — > u n+ \ = u n ± A and u n+ i — > u n+ 2 = Un+i ± A thus getting u n — > u n+ 2 = 
Un ± 2 • A, together with two times u n — > u n+ 2 — u n . Considering the coefficients of the biquadratic 
( ff7|) one thus gets, very simply, a polynomial representation of the shift doubling A — > 2 • A : 

7 > 7"( 2 ) — _ 7 2 7 2 1 7 2 7 2 _ 7 2 7 2 

J i J x J z J i T z °y x °y 

7 > 7( 2 ) - 7 2 7 2 _ 7 2 7 2 _ 7 2 7 2 

Jy ,J y — J ,J X J ,J y J x ,J y 

J Z ► ' = —J z J x —J z Jy J T Jr,. Jy (49) 



The modulus (46) is (as it should) invariant by (ft9j), which represents the shift doubling transformation. 



Of course there is nothing specific with the shift doubling : similar calculations can be performed to 
get polynomial representations of A — > M ■ A, for any integer M . A general calculation, corresponding 
to eliminations between two biquadratics Ti of same modulus (|46|), but different shifts A and A', will 
be given elsewhere. 

Terminology problem : From the point of view of discrete dynamical systems a mapping like ( f49| ) 
(and the mappings given below (Ell), (|54|), ...) could, at first sight, be called "integrable" : the iteration 
of this (two-dimensional) mapping "densities" algebraic curves (elliptic curves) foliatin g th e whole two- 
dimensional space, namely M d = constant, exactly as an integrable mapping does [|l4[|l6f |. One can 
even write explicit analytical expressions for the iV-th iterate, for any N . However, this mapping is 
not reversible, the growth of the calculations [pl| , p2| , ^3| is exponential ( 2 N exponential growth, ln(2) 
topological entropy, ...). In fact this very example of "calculable" chaos is the exact equivalent of the 
situation encountered with the logistic map x — > a • x • (1 — x) for a = 4 : one does not have a 
representation of a translation 9 — > 9 + N ■ A, but a representation of the iteration of a multiplication 
by 2 : 9 -> 2 N ■ 9. 



VI. POLYNOMIAL REPRESENTATIONS OF THE MULTIPLICATION OF THE SHIFT BY 
AN INTEGER AND ASSOCIATED FINITE ORDER CONDITIONS 



The multiplication of the shift by three can be obtained using the previous elimination procedure, 
namely eliminating y between ^(x, y) and T\{y, z) (or equivalently eliminating y between Ti(x, y) 
and 1^2(2/, z)), thus yielding a resultant which factorizes into two biquadratics of the same form as the 
two previous ones, namely T\(x, z) and T^(x, z) : 

T 3 (x, z) = (x 2 - 1) ■ (z 2 - 1) • 4 3 ) - (x 2 + 1) ■ (z 2 + 1) ■ Jf ) + 4 • Jj 3 ) • x ■ z = (50) 

where ji 3 \ Jy and J^ are given below. 

The multiplication of the shift by three has the following polynomial representation (J x , J y , J z ) — > 

I t(3) t(3) 7 (3)s . 
(Jx , Jy , Jz ) ■ 

J X = j X ■ (— 2 j z jy J x — 3 jy j Z + 2 jy ^ ^ + ^ ^ ^ ^ ^ ^ ^ + ^ ^ ) 

= jy ■ j;, J X + 2 j;, jy J x j rjyj x — 2 jy ^ ^ + ^ ^ +2 ^ ^ ) 

Jf) = J z ■ ( J a 4 J, 4 + 2 J 4 J z 2 J 2 - 3 Jl Jt - 2 J 2 Jt J 2 + 2 J 2 J 2 Jt + Jt Jt) (51) 
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The modulus (Q) is (as it should) invariant by the polynomial representation (|5l]) of the multiplication 
of the shift by three (J5l]) . 

The multiplication of the shift by four has the following polynomial representation (J x , J y , J z ) — > 
(J (4) J (4) J (4) ) • 

J x -* = —4 jj J y J z — 6 J z J y J x + 4 J z J y Jj, + 4 J 2 J y jj + 4 J z J y J x — A J z J x J y + 10 J z J x J y 

— 4 — j^, jy — J Z jy — j Z j^. + 4 jy J x — A j z J x jy + A ^ ^ ^ ~ 6 ^J. ^ 

^ = —j x jy—j z jy—j z j X —A j jy J z — 6 j g j y j Q. ^ ^ ^J, ^ ^J. ^" 4 ^ ~ ^ ^ ^ 

+ 10 jy j^, — A j z jy j^, — A j z jy jy. A ^ ^ ^ ~ 6 ^ ^J. ^ ^" 4 ^ ^J, ^ ^ ^ ^ ^ ^ 

J z ^ = —J x jy—J z jy—J z j x ~6 j z J x jy + A j g j rj , j y ^ A j g jy j 3 ~ 4 ^ ^ 4" 10 ^J, jy j^j, 

-A J* Jl Jl - A Jl J 2 y J% + A Jl J e y J 4 - A Jl Jl Jy + A J 2 Jl Jl + A J\ J* Jl - 6 J 2 4 J 8 y J 4 (52) 



which can be obtained, either by the elimination of y between ^(x, y) and ^(y, z) (and extracting a 
{x— z) 2 factor in the resultant), or, equivalently, by the elimination of y between T\(x, y) and T^(y, z), 
or the elimination of y between ^(x, y) and T\(y, z) (and extracting a T2 factor in the resultant). 
Again, one gets T^(x, z) : 

T 4 (x,z) = (x 2 -l)-{z 2 -l)-J^-(x 2 + l)-{z 2 + 1)-J^ + 4- jW-x-z = (53) 



where Jj , Jp and J z 4 ^ are given above. It can easily be verified that ( |52] ) can be obtained directly 
combining (B9[) with itself. 

The multiplication of the shift by five has the following polynomial representation {J x , J yi J z ) —> 

( t(5) t(5) x(5)s . 

Jx ^ — -^x ^(.Jxi Jyi Jz) 

Jy ^ — Jy ' Py ^{Jxi Jyi J z) ~ Jy ' Pg° ^(.Jyi Jz"i Jx) 

j| 5) = J z ■ lf\j x , Jy, Jz) = Jz ■ ^i 5) (J Z , J*, Jy) (54) 



where : 

Ph[Jx\ Jy, Jz) — J z J x ^ ^ Jz Jy + 6 Jy J z J x + 14 J y J z J x — 29 J y J z J x — 6 J y J z J x 

-20 Jl Jl Jf + 15 Jl J 4 J x 2 - 36 Jl Jf Jl + 50 Jl Jf J* - 9 jj J 12 J 4 + 15 J 4 jj J^ 2 

— 20 J Z J X +4 jy jj, jj. + 36 jy j^, j^. ~ 20 jy j^, j^. ~ 10 j y j ^ j^, 4" 34 jy J z jj, 
+ 14 jy j z j x +4 jy j z j x 6 jy j z j^j. 36 jy j z j^. 

— 9 jy jj, jj. — 29 jy j^ jj. — 10 jy jj, j^j, + j y j^. ^ ^ j y j % j ^ j y ^ ^ j^ (^5) 

The modulus ( |4^ ) is, again, invariant by the polynomial representation of the multiplication of the 
shift by five, namely (|54|). One remarks that Px{J x , Jy, Jz) singles out J x and is invariant under 
the permutation J y <-> J z and, similarly, Py 5 \j Xl J yi J z ) singles out J y and is invariant under the 

permutation J x «-> J z and Pz 5 \j x , J y , J z ) singles out J z and is invariant under the permutation 
J x <-> Jy. One has similar results for the polynomial representation of the multiplication of the shift by 
M = 6, 7, 9, 11, • • • The explicit expressions of these polynomial representations are given in Appendix 
A. 

Let us denote by Tjv a biquadratic corresponding to u — > u ± N ■ A : 
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T N (x,z) = {x 2 -l)-(z 2 - 1) • jW - (z 2 + l)-(z 2 + !)• Jf»+ 4- jM.j.z = (56) 



In general, it should be noticed that the elimination of y between Tm{x, y) and Tm'{v, z), yields a 
resultant which is factorized into TrM+M')( x > z) and TrM-M')( x > z ) (f° r M > M'). When seeking for a 
new Tn(x, z) there may be many (M, M') enabling to get Tn{%, z) (that is such that N = M + M'). 
One can verify that all these calculations give, as it should, the same result (in agreement with a 
polynomial representation of u — > u ± M ■ A ± A/' • A giving u ± (M + M') - A or : u ± (M — M') • A. 
Let us denote TV these homogeneous polynomial representations of the multiplication of the shift by 
the natural integer N. In the same spirit one can verify, for N = M ■ M' (N, M, M' natural integers), 
that : 

Tpf(J x , Jy, J z ) — (Tm) M (Jx, Jy, Jz) = Tm(Tm(Tm(- ■ -Tm(Jx, Jy, Jz)'"))) = 

= (T M >) M (Jx, J y , Jz) = T m ,(T m >(Tm>(- ■ ■ T M >(J X , J y , J z ) ■ ■ ■))) (57) 

One can, for instance, easily verify that T2 and T3 commute, as well as T2 and T5. Similarly one can 
verify, in a brute-force way, that T3 and T5 commute. This commutation relations are true for TV and 
Tjv/, for any N and M. One thus has a polynomial representation of the natural integers together with 
their multiplication. 

In general, the multiplication of the shift by a prime number N ^ 2 has the following polynomial 
representation (J x , J y , J z ) — > (Jx N \ Jy N \ Jz N%> ) ■ 

J ■ p( n ) ( J J T ) 
• J y 1 x \ ,J yj ,J z, °x) 

Jz-P { x N) {Jz, Jx, Jy) (58) 

One remarks that the Px"\Jx, Jy, T z )'s (and thus Py"\j x , J y , J z ), Pz N \j x , Jy, Jz)) are polynomials 
of J 2 , J 2 and J 2 . One verifies easily that the homogeneous polynomial transformations Tm are all of 
degree M 2 , in J x , J y , J z , for M = 2, 3, 4, 5, 6, 7, 9, 11. 



j{N) _ j . p(N)(j j j \ 

°x ,J x 1 x \' J x, -Jy, •Jz) 

Jy ' = Jy' Py \jx, Jy, Jz) = 

J z ^ — Jz' P z ^{Jx, Jy, J z) — 



A. Finite order conditions and associated algebraic varieties 



Let us show that one can deduce the (projective) finite order conditions K M (R) — ( ■ R, from the 
previous polynomial representations. 

Since one knows that the (projecti ve) fin ite order conditions of K 2 often play a singled-out role for 
integrability (see the previous section (ElIAD ), and, in particular, since one knows B that the free-fermion 



conditions of the Baxter model correspond to K A (R) = (^-R, one can, as an exercise, try to systematically 
write, for the sixteen vertex model, the (projective) finite order conditions K 2N (R) = £ • R, with N 
natural integer, the corresponding algebraic varieties being "good candidates" for new free-(para?)- 
fermions, or new equivalent of the integrable chiral Potts model [0,^O). Let us first give these finite order 
conditions for the Baxter model where these algebraic varieties are naturally associated with the set of 
RSOS models @. 

Let us recall @j that the free-fermion condition for the XYZ Hamiltonian, J z = 0, corresponds to the 
finite order (projective) condition K 4 (R) = £ • R. Recalling the polynomial representation (j4TJ) of the 
shift doubling, one can easily get convinced that J^ — should correspond 111 to K 8 (R) — £-i?. This 



12 The idea here is the following : K, or K, corresponding, with some well-suited spectral parameter, to 
8 — ► 6 + 77, K 2 , or K 2 , must correspond to 9 — > 8 + 2 ■ 77 . A finite order condition of order M 
corresponds to a commensuration of n with a period of the elliptic curves : n = V/M. Changing K into K 2 
amounts to changing n into 2 • n, or equivalently, changing the order M into 2 M. 
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can be verified by a straight calculation. More generally, a polynomial condition Cm(Jx, Jy, J z ) = 0, 
corresponding to K (R) = ( ■ R, has to be compatible with the polynomial representations of A — ► 
N ■ A. This compatibility is, in fact, an efficient way to get these finite order conditions. Explicit 
expressions of finite order conditions, as well as their compatibility with the shift doubling and, more 
generally, the polynomial representations of A — ► M ■ A, are given in Appendix B. For a prime integer 

N ^ 2 the algebraic varieties P^ N) (J X , J y , J z ) = 0, P^ N) {J X , J v , J z ) = 0, and P { z N \j Xl J v , J z ) = 
give order 4 N conditions : 

K 4N (R) = (-R (59) 



Since the P^'s (and the J^'s for N even) are functions of Jj, J 2 and J 2 , the order AN conditions, 
K 4N (R) = X ■ R, are also functions of the square J|, J 2 and J 2 . One can easily get infinite families 
of finite order conditions. For instance, iterating the shift doubling ( pt9| ) (resp. (|5l|)), and using this 
transformation on J z = 0, one easily gets an infinite number of algebraic varieties corresponding to the 
finite order conditions of order 2 N (resp. 3^). Combining (^9|) and (|5l]), one gets straightforwardly the 
finite order conditions of order 2 N x 3 M . More details will be given elsewhere. 

Remark : Finite order conditions for the sixteen vertex model. 

Using the previous results, the decomposition ( |39"|) of the sixteen vertex model into a "-fT 2 -effective" 
Baxter model, and relation (f4(i|), one can obtain the finite order conditions, K AN (R) = £ • R, for 
the sixteen vertex model and find that they are actually codimension-one algebraic varieties. Recalling 
( p0| ) , one easily deduces that these finite order conditions are actually given in terms of the finite order 
conditions of the "if 2 -effective" Baxter model. These finite order conditions are simply expressed in 
terms of the associated K 2 -effective" variables J x , J y and J z which can be obtained from relations 
©, ©■ 

From this example one sees that the, at first sight, "hardly Baxterisable" case of finite order iteration 
provides, to some extend, more results, and structures, than a "standard" infinite order Baxterisation: 
one gets, for instance, a polynomial representation of the natural integers together with their multiplica- 
tion, this polynomial representation leaving invariant the modular invariant p6| of the elliptic curves, 
and giving codimension-one algebraic varieties compatible with this structure ... 



VII. LET US BAXTERISE QUANTUM HAMILTONIANS 



Let us now consider a typical problem for any theoretician who wants to provide some "interesting 
contribution" in High-Tc. Let us consider a str ong ly correlated quantum Hamiltonian which looks like 
a t-J model |37|j38|j3g[| ) or a Hubbard model jl(^,[ll||, or some coupled XYZ quantum chains : how is it 
possible to see if this quantum Hamiltonian is integrable ? More generally, let us consider a quantum 
Hamiltonian. How to see if it is possible to solve this quantum Hamiltonian ? The Bethe Ansatz only 
works if one has some "conservation operator somewhere" L3 that enables to see this model as some 
"avatar" of the six- vertex model (XXZ chain) . Beyond this restricted "six- vertex" framework, and if one 
cannot, or does not know how to, associate a (commuting) family of transfer matrices, commuting with 
this Hamiltonian, one has very few tools left to solve this quantum Hamiltonian Ej. For a Yang-Baxter 
integrable model depending on one spectral parameter, the associated integrable quantum Hamiltonian 
can be seen as the derivative of the (logarithm) of the transfer matrix at a singled-out value of the 
spectral parameter : therefore the integrable quantum Hamiltonian does not depend on the spectral 
parameter. If the Bethe Ansatz is too complicated, or does not exist (higher genus curves, ...), how can 



13 The typical example is the operator '^2 n (J n f° r the XXZ Hamiltonian p^ , [t3{ . 

14 The random matrix theory analysis of the level spacing distribution remains a possible tool to simply "de- 
tect" integrability J44|]. However, besides the technical difficulties associated with the unfolding procedure, the 
calculations become very large for the 16 x 16 (Hubbard) i?-matrices, or for two coupled spin chains. 
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we make the spectral parameter (s) "emerge" so that the integrability structure becomes crystal clear ? 
This is clearly another type of Baxterisation problem. Let us sketch how this can be done. 



A. Let us Baxterise the XYZ quantum Hamiltonian 



Let us recall the XYZ Hamiltonian p9| , p5[ |, Hxyz = ^ n H n ,n+i, and let us represent H n ^ n+ i = 
Jx • &n • &n+i + Jy ' a n ' a n+i + ^ z ' G n ' a n+i ano - ^ i the matrix of permutation of the vertical space n 
and n + 1 , as 4 x 4 matrices : 

Jz 



H 



71,71+1 



Jx- Jy 

-Jz Jx + Jy 

Jx + Jy -Jz 

Jx- Jy Jz 



P 



10 
10 
10 
1 



(60) 



Recalling Baxter's notations, the 7?-matrix, reads at order one in some small expansion parameter e 



R= P + e-P-H n 



n+l 



1 + eJz 



e- (Jx - Jy) 



- 1 








0" 










1 










1 










. 








1 . 







e • ( Jx + Jy) 
1 - e • Jz 




1 





e • Jz 



e ( Jx + Jy) 




Jz 



Jx — Jy 
(Jx - Jy) 



1 + e ■ Jz 



Jx - Jy 

Jx + Jy —Jz 
— Jz Jx + Jy 
Jz 

a d 
b c 
c b 
d a 



(61) 



One has the correspondence : 
a = 1 + e ■ J z + ■ ■ ■ , b = e ■ (J x + J y ) + 



l-e-J z 



e ■ ( J x 



One verifies, immediately, that the algebraic i^-covariants of the Baxter model read : 
a 2 + b 2 - c 2 - d 2 

= 2e • J z , a ■ b + c ■ d — 2e ■ J x , a ■ b — c ■ d — 2e ■ J y 



(62) 



the if-invariants being (for instance) the ratio J x /J z and J y /J z . The i?-matrix P can be seen to belong 
to all the elliptic curves of this foliation of the Baxter model (the invariants (a 2 + b 2 — c 2 — d 2 )/a/b 
or ab/c/d are of the form 0/0). Matrix P is really the equivalent of the base point of an elliptic 
foliation @. 

Having the XYZ Hamiltonian and willing to "recover" the Baxter i?-matrix (and its canonical el- 
liptic parameterization) is a slight modification of the previous Baxterisation process of i?-matrices, 
or monodromy matrices, where we were building the algebraic variety from one arbitrary point of the 
algebraic variety. We do not have a point of the elliptic curve here, but rather a singular point, which 
is the equivalent of the base point of an elliptic foliation, namely point P, and a "vector" ( J x , J y , J z ) 
giving the tangent to the elliptic curve at point P. It is however clear that [P, (J x , J y , J z )\ is sufficient 
to build, in a unique way, the elliptic curve, and thus the spectral parameter. 

A heuristic remark : Point P is a singular point of K 2 and the diagonal matrices are fixed points 
of K 2 : 
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K(P) 
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More generally, considering the most general 4x4 matrices (sixteen vertex model), the singular matrices 
for K 2 are the i?-matrices such that K{R) is a rank two matrix. For the Baxter model singular matrices, 
such that K 2 (R) is the null matrix, are, for instance, d = b = 0ora = c = 0. The base points of the 
elliptic foliation of the parameter space of the Baxter model are a subset of these singular subvarieties. 
The base points of the elliptic foliation are the points such that the if- invariants, (a 2 + b 2 — c 2 — d 2 ) / a/b 
or a b/c/d, are both of the form 0/0, namely d = b = 0, c = a, or d = b = 0, a = — c, or a = d = 0, b = c, 
or b = — c, a = d = 0. This can be easily ge nerali zed to the sixteen vertex model by imposing that the 
PiS if-covariants Pq] , mentioned in section ( V A ), are all equal to zero. 



In general these two sets of subvarieties, namely the singular subvarieties and the subvarieties of fixed 
points of K 2 , play a crucial role in the Baxterisation of a Hamiltonian. The subvariety of fixed points 
of K 2 always contains the set of all the diagonal matrices. Not surprisingly we will see, in the next 
section, that diagonal matrices naturally occur (in a non-trivial way) in the Baxterisation process (see 
for instance (69) below). 



B. Let us Baxterise the t — J quantum Hamiltonian 



Let us recall the Hamiltonian of the t — J model |37| , and let us consider an ordering of the 9x9 
i?-matrices, well-suited for the partial transposition t\ : 



(+,+), (+,-), (+,0), (-,+), (-,-), (-,0), (0,+), (0,-), (0,0) 



(63) 



With this ordering, the "local Hamiltonian" of the t — J model, equivalent of the previous H n>n +\ , and 
the permutation matrix P, read respectively : 
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The integrability cases (47| correspond to (V, J) = (-t/2,2t) or (t/2, -2t) or (3t/2,2t) or 
(— 3i/2, — 2t). The i?-matrix reads at order one in some small expansion parameter e: 



R = P + e-P- H n , n+1 
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The (supersymmetric) integrable case (V, J) 
tonian" H nn +i( J = 2 1, V = - t/2) 



= (—t/2, 2t) yields the (supersymmetric) "localHam.il- 



t ' H susy t 



t 


H S usy • 












































-1 





1 



































-1 











1 





-1 

































































-1 











-1 



































-1 







































(66) 



Let us "absorb" the homogeneity of H susy by introducing x 
Baxterise: 



Ra,„ = P-X-P-H, 



susy 



+ 



Let us introduce the two 9x9 matrices 



N = 



-e ■ t. At first sight, one wants to 

(67) 



N' = 



- 1 
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■ 1 
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(68) 



With ordering (63) matrix N is associated with the electron- counting operator which commutes with the 



t-J Hamiltonian. Matrix TV commutes with P and is such that the matrix Di ag — P-H susy — P + P-N 
is a diagonal matrix of successive diagonal entries : 1, 1, —1, 1, 1, —1, —1, —1, —1. Actually, as far as 
the integrability of the model is concerned, one can add to the Hamiltonian an operator which commutes 
with it. Obviously one can also add to the i?-matrix, a term x ■ P which just amount to a global I + x 
multiplicative factor in front of P. This amounts to introducing, instead of (p7]), the following i?-matrix: 
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Ri(x) = P - x -D. 



lag 



(1 



{l+x)-P- (id 



■P- x-P-H 



susy 



(H susy +N)\ = 



- X ■ 
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1 + x 



(69) 



where Id denotes the 9x9 identity matrix. Under transformation the R- matrix Ri (x) becomes 
another matrix, of the same form Ed, but where x — > x + 2 : 



K A (Ri{x)) 



1 



Ri(x + 2) 



(70) 



This solves the Baxterisation problem of (|66|), reducing it, after addition of "well-suited" commuting 
operators, to a simple linear interpolation. 



• Another integrable case J47j for (64) is (V, J) = {it/2, 2t), for which one has to introduce another 

the 



for 



correction namely N' in (p8[), yielding the same i?-matrix (|6Sj) , namely R\ = P — x ■ D 
diagonal matrix Di ag being, now, equal to P ■ H t .j — P + P ■ N', where H t j is given by 
(V, J)= (3t/2,2t). 

• Two other integrable cases @ for @ namely (V, J) = (-3f/2, -2t) and (V, J) = (t/2, -2t) 
yield the same i?-matrix i?2- For (V, J) — (—3t/2, —2t) and (V, J) = (t/2, —2t) one has to introduce 
the correction N and N' respectively, yielding a very simple i?-matrix R2 = x-Id + P- Similarly, under 
transformation K 2 , matrix R2(x) becomes another matrix of the same form but where x —* x — 3 : 



K 2 (R 2 {x)) 



1 



x ■ (x — 3) 



iJa(s-3) 



(71) 



thus solving the Baxterisation problem of these last two integrable cases as a simple linear interpolation 
between Id and P. 



C. Hamiltonian dependence of the Baxterisation procedure 



One certainly wants the Baxterisation procedure of a quantum Hamiltonian to be "universal" : it 
should be compatible with the addition of quantum operators which commute with the quantum Hamil- 
tonian, and, thus, do not modify the integrability of the Hamiltonian. 

The Baxterisation procedure is, of course, compatible with the gauge transformations on the quantum 
Hamiltonian : 



H, 



n,n+l 



(72) 



15 Note that K 2 (Ri) has a form similar to i?i(a;-|-l) up to some change of signs of some entries (supersymmetric 
graduation) . 
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where X n and I n +x denote the identity operators on site n and n+1 respectively, the "exponentiation" 
of (72) giving the R — > ® h^ 1 ■ R ■ g ® h well-known gauge symmetries [K4|, compatible with the 



Baxterisation: K {g ®h Rg <8> h) = <7 _1 <g>/i _1 -K 2 (R)-g ® ft,. Furthermore, from the symmetries of 
the Baxterisation (see for instance (|32|)), one deduces that the Baxterisation of a Quantum Hamiltonian 
is necessarily compatible with the following deformations of H n ^ n+ i generalizing E3 ( ^2[) : 

-> #r»,n+l + (G" ® X„ - Z n+1 <g> G£ +1 ) + (G£ ® X„ + T n+1 ® G^ +1 ) (73) 

The Baxterisation should also provide (at least as far as integrability is concerned) compatible results 
when one adds, to the quantum Hamiltonian ^2 n H nin +i, operators commuting with it : 

H n , n +i — * H n<n+ x + 0" n+1 , a= l,---,r (74) 

where the J2n^n n+i' s commute with J2n ^n,n+i- When the quantum Hamiltonian J2 n H n ,n+i is 
integrable one would like to have a complete description of all the integrable deformations (fri)), and a 
complete description of the algebraic variety associated with the Baxterisation of the largest integrable 
deformation of the quantum Hamiltonian (as many O™ „ +1 's as possible). Are all the integrable de- 
formations O" _ , 1 necessar ily of the ^-compatible form @ ? All the 0° n+1 's are such that the 
S n G" n +i's commute with ^ n ff n , n +i, but do all the O" n+1 's commute with H, ltn+ i and commute til 
all together ? Is there a way to find r, the largest number jif integrable deformations O™ „ +1 , by simple 
arguments? Is it possible to find lower and upper bounds Il3 for r ? All these questions will be analysed 
in details elsewhere. Let us just give, in the following sections, some partial answers to these questions 
based on deformations of the two previous "heuristic" integrable cases, namely i?-matrices Ri(x) and 
R 2 (x). 

Let us try to see what happens when one changes the integrable Hamiltonian adding an operator 
which commutes with it. Suppose that, instead of introducing the i?-matrix Ri(x), which yields a 
simple linear interpolation for the Baxterisation process (|69[), one introduces : 

R(x, y) = P - x ■ D lag + y-P-N = (1 + x) ■ P - x ■ P ■ H susy - (x - y) ■ P ■ N + ■ ■ ■ (75) 

One can actually easily verify that the n-th if 2 -iterate of R(x, y) are all linear combinaisons of R(x, y), 
K 2 {R{x, y)), K\R{x, y)), K e (R(x, y)) and K 8 {R{x, y)) or, equivalently, R{x, y), K 2 {R{x, y)), 
K 4 (R(x, y)), P and P N : 

K 2n (R(x, y)) = (76) 

= A { n) ■ R(x, y) + A[ n) ■ K 2 (R(x, y)) + A { 2 n) ■ K\R(x, y)) + A™ ■ P + A { 4 n) ■ P ■ N = 

= B { n) ■ R(x, y) + B[ n) ■ K 2 (R(x, y)) + B« ■ K\R(x, y)) + flW • K 6 (R(x, y)) + B™ ■ K 8 (R(x, y)) 

In other words, the Baxterisation acts in the four-dimensional vector space of the i?-matrices of the 
form ([76]). The Baxterisation of R(x, y) gives successive 9x9 matrices K 2n (R(x, y)) of the form: 



16 A simple example corresponds, for instance, to adding operator &n to the XX Z Hamiltonian [ j^2|j43| . 

17 In the case where K 2 is an infinite order transformation densifying an algebraic variety, this variety is an 
Abelian variety and one can deduce, from this Abelian property, that the P ■ 0™ !jl+1 's should commute. 

18 For Abelian varieties parameterized by theta functions of g variables, one could imagine that r should be 
related to g + N g , where N g denotes the number of independent gauge deformations (fr2[). 

19 Five homogeneous parameters. Note that the odd jf-iterate do not belong to this four-dimensional vector 
space but to its transformed under ti. 
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R(A, B, C, D, E, F) 
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where : A = B + C 



(77) 



Factorization scheme : Let us Baxterise too = R(x, y) using the homogeneous transformation K. 
The first iteration reads toi = K(mo). Let us denote G\ the gcd of all the entries of matrix m\. In 
order to go on iterating, it is better to introduce the "reduced" matrix, Mi — m\/G\. Let us denote 
to 2 = K(M\). Similarly, one can introduce G 2 the gcd of all the entries of matrix m%, and define the 
"reduced" matrix, M 2 = m 2 /G 2 , and so on. Introducing G(w) the generating function of the degree 
of the successive gcd's G n 's, and D(w) the generating function of the degree of the determinants of the 
successive reduced matrices M ra 's, one gets for the iteration of (for instance) R(x, 2) : 

G(w) = 4w + 25w; 2 + 34u; 3 + 144u; 4 + 164u; 5 + 604w; 6 + 684 w 7 + +2444w 8 + 2764w 9 --- (78) 
D(w) = 9 + 36w + 63w 2 + 198w 3 + 288u; 4 + 828w 5 + 1188w 6 + 3348w 7 +4788w 8 + 13428w 9 + ••• 



Of course one easily verifies that (1 — 8w) ■ D(w) + 9G(w) — 9=0, which is simply deduced from 

K(M n ) 



A I. 



n+l 



G n 



+i 



(79) 



One easily verifies, up to order nine, that ( |78| ) are the expansions of the two rational expressions : 
w ■ (4 + 21 w - 7 w 2 +26w 3 - 16 w 4 ) 



G(w) 



(1 -w)(l + 2w) (1 -2w) 



y ' (l-w)(l + 2w) (1 -2w) K ' 



There is nothing specific with y = 2 : beyond R(x, 2), one also gets (generically) a 2™ exponential 
growth of the iteration calculations for R(x. y), thus excluding Yang-Baxter integrability. This may, 
at first sight, seem in contradiction with the fact that the (local) Hamiltonian H = H susy + p ■ N 



is "as integrable as" H susy (or H, 



L susy 



N). In fact the (local) Hamiltonian H = H, 



s usy 



■ N 

should also yield a polynomial growth of the calculation corresponding to some algebraic subvariety of 
((77]) such that the tangent space to this algebraic subvariety, at point R = P. contains the vector 
P ■ H = P ■ (H susy + p ■ N). In other words, R(x, y) should only be integrable in the (x , y) — > (0, 0) 
limit, but not for finite values for [x , y). In our case the plane containing point P and the two integrable 
"vectors" P ■ H susy and P ■ (H susy + p ■ N) , is not included in the algebraic subvariety of polynomial 
growth. This is a quite general situation : in general, integrability does not correspond to linear spaces 
but to algebraic subvarieties with some "curvature" . 



Deforming ( |69| ) with P ■ N, one thus finds that the integrability of (|69|) is "surrounded" by a (four- 
dimensional) space corresponding, generically, to a 2™ exponential growth. However this does not mean 
that the whole four-dimensional space (^) (except ( |69|) o f course) corresponds to a 2™ exponential 



growth: it is possible that some algebraic subvarieties of lYTA) could correspond to a polynomial growth of 
the iteration calculations and, possibly, correspond to a Yang-Baxter integrability. The question is thus 
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to find such polynomial growth, and possibly Yang-Baxter integrable, algebraic subvarieties l3. The next 
section gives an explicit example of particular Yang-Baxter integrable deformations of the simple linear 
interpolation R2 = Id + x-P, and gives the associated integrable algebraic subvariety. Section (VIIE) 
sketches a strategy to find systematically Ej these polynomial growth, and integrable, subvarieties (or, 
more generally, the subvarieties of smaller topological entropy pi 22 1). 



D. Beyond simple linear interpolation : an integrable deformation of R2 — Id + x ■ P 



Matrix R2 = Id + x ■ P gives a Baxterisation of some integrable t-J Hamiltonian as a simple linear 
interpolation between P and the identity matrix Id- Let us try to extend this simple linear interpolation, 
and "merge" it into a larger integrable i?-matrix. Similarly to what has been previously performed, let 
us consider two matrices N ew and A sym , which respectively commute, and anticommute, with P : 
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Do note that P ■ A sym is a diagonal matrix with diagonal entries 0, —1, 1, 1, 0, —1, —1, 1, 0. Let us 
consider a deformation of R2 = Id + x ■ P with these two matrices. This amounts to considering 
R = P + a ■ Id + (3 ■ A sym + 7 • N ew , or equivalently the non-symmetric i?-matrix : 



R z (t, z) 



This form (g3|) is stable by h 
well-suited for the action of K 2 . Let us consider 
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(83) 



algebra is thus an algebra 



for : 



This should be particularly difficult when these subvarieties are not codimension-one subvarieties, but codi- 
mension two, three (or more ...) subvarieties 

21 Of course one can always iterate numerically a point, very close to point R = P, moving away from this 
point along directions like (for Ri) P ■ H SU s y , P ■ N or P ■ N' and visualize these integrable subvarieties [ [l2| . 
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A = t • (1 - z 6 ) , B = z 4 ■ (t 2 - 1) , C = z 2 ■ (t 2 - 1) , D = t 2 — z 6 (84) 
The corresponding i?-matrix (§3)) Baxterises in a very simple way : 

f3(t, z) ■ K 2 (R z {t, zf\ = R z {t,t-z) (85) 

(1 - z 3 ) (1 + z 3 ) (zH 3 - 1) (z 3 t 3 + 1) t 2 

where : 0(t, z) = V u > \, , \, , ^ ^7- ^f— 

y ' (t + z 3 ) (z 3 t + 1) (z 3 t - 1) (t - z 3 ) 

The i?-matrix (]83|), with (|84|), is actually one of the Yang-Baxter integrablc models of Perk and 
Schultz and can be seen as a deformation of Ri(x). Actually, in the limit (z, t) = (1 + 

Z ■ h, 1 + T ■ h), where h — > 0, one finds that operator A sym occurs at order one in h , while operator 
N ew occurs at order two in h : 



z 2 -(z 2 + l)-{t 2 ~l) \ T J \ " y "' 2T 

f Z-(8Z-3T) Z 2 Z(16Z 2 -3TZ-3T 2 ) x 2 

-I- ^ 2 ' ew 2~ ' sym 47* d ) 



Parametrization (|84|) is equivalent to the (Yang-Baxter) integrability j^,|49j condition : 

C ond (A, B, C, D) = B 3 D - B 2 C 2 + C 3 D + BCA 2 - BCD 2 = (86) 

This algebraic variety is an extension of the simple interpolation between the identity matrix and P 
(Yang R- matrix), which corresponds to C = B and D = A + B. Actually, for C = B, condition 
(|S6|) factorizes as follows : 



C nd(A, B, B, D) = —B ■ [D — (B — A)j ■ [D — (A + B)j = (87) 

Considering the general i?-matrix of the form j83| ) (beyond the integrability condition ((8^)), the 
factorization analysis, which extracts the gcd's at every step (see (^9|)), yields a 2" exponential growth 
of the calculations. The generating functions of the degrees of the M n 's and G„'s (see (|79|), (|78|)) read 
(the expansions of D(w) and G(w) have been obtained up to w 8 ) : 

s „ 3-2w + 6w 2 + 2w 3 -6w 4 „. . 6 - 7w + 18w 2 + 6w 3 - 16 w 4 

(1 — w) (1 + w) (1 — 2 w) (1 — w) (1 + w) (1 — 2w) 

In order to have a polynomial growth of the calculations some additional factorizations must occur. 
Actually, fop-([B4|), that is when the integrability condition (|8^) is satisfied, the generating function of 
the degrees E3 reads (the expansions of d(w) and g(w) have been obtained up to w & ) : 

n 6-2w + 7w 2 + 7w 3 . . (l4-9w + 19w 2 + 18w 3 ) 

d[w) = 3 — r — , g(w) = w j— r 89 

(1 — w z ) (1 — w) (1 — vr) (1 — w) 



22 The degree, here, being the degree in t, the calculations being performed for ( p4| ) with z = 3. Note that 
there is nothing specific with z — 3. 
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E. Integrability emerging from the occurrence of additional factorizations 



In fact one can actually find the integrability condition (|36|), as a condition corresponding to the 
occurrence of additional factorizations. The "bifurcation" from the polynomial growth of the degrees to 
a 2™ exponential growth of the degrees, takes place immediately when leaving the integrability condition 
(|S6|). In the (A, B, C, D) parameter space, the integrability condition ( p^ ) is "surrounded" by a 2" 
exponential growth. More specifically, the "bifurcation" from the polynomial growth of the degrees to 
a 2™ exponential growth of the degrees, takes place with the first gcd, namely G\. Let us consider 
the factorization of the i?-matrix corresponding to (84) (that is (|86|)), and the one where A has been 
replaced by A + u, which corresponds to a 2™ exponential growth. The first gcd's corresponding to these 
two situations, namely for u = Q_ (polynomial growth, integrable) g\, and for u ^ (2 n exponential 
growth) G\(u), read respectively : 

Gx(u) = (-u + z 3 - t + tz 6 - zh 2 ) 2 -(-u-z 3 -t + tz 6 + zH 2 ) 2 ■ (t 2 - z 6 ) 2 = D 2 -(B-C - A 2 ) 2 

9i = (t 2 -z 6 ) 5 -(zh 2 -!) 2 (90) 

One remarks that : 



9i 



(i 2 -z 6 )-Gi(u = 0) 



D ■ G x {u = 0) 



(91) 



which amounts to saying that, restricting to the integrability condition ( |86| ) (here u = 0), an additional 
factorization occurs, changing (jS8|), and its 2 N exponential growth of the calculations, into (|89|) and 
its associated polynomial growth. When the integrability condition ( |S6| ) is not satisfied, one has G\ = 
D 2 ■ (B ■ C — A 2 ) 2 , and the first "reduced" matrix Mi reads : 







'A 2 


-BC 











DC 











-DB ' 








DA 



































DA 






















K(R) 













DA 

















Mi = 




DB 








A 2 


— BC 











-DC 




Gi 











































DA 



































DA 



































DA 











DC 











DB 











A 2 — BC _ 


Condition (| 


(]), namely (B 3 - 


h C 3 )D 


+ ( 


'{A 2 


- BC) - 


D^BC = 


0, is 


nothing but the 



(92) 



divisibility of A 2 — BC by D, which obviously enables to extract an additional factor D in all the entries 
of (9^). The fact that (86) is such a condition of divisibility, becomes obvious if one uses the fact that 



(36) is a rational surface, and one actually uses this rational parameterization. With parameterization 
(~) the division of A 2 - BC by D reads : 



BC 



D 



= 1 



z 6 -t 2 



(93) 



23 Note that the factors of G\(u) actually correspond to some factors in &et{R(u)) which reads &et(R(u)) = 
(-u + z 3 - t + tz 6 - z 3 t 2 ) 3 ■ (-u - z 3 - t + tz 6 + z a t 2 ) 3 ■ (t 2 - z 6 ) 3 . This is a general result : the cofactors 
associated with K" or ti and I are necessarily related with powers of the factors occurring in the determinant 
(or the whole determinant), and their transformed by the group generated by ti and I. 
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Finding systematically algebraic varieties (like ( pq ) but not necessarily rational) which correspond 
to additional factorizations in the factorization scheme, and thus correspond to a smaller complexity 
of the calculations (smaller Arnold complexity , smaller topological entropy p2] | ) , can be seen as a 
new approach of integrable models. The set of all these singled-out algebraic varieties, of non generic 
topological entropy pi] ], should be a complicated stratified space. A remarkable subset is the set of 
algebraic varieties associated with a polynomial growth of the calculations. A remarkable subset of the 
previous subset is the set of algebraic varieties which are Yang-Baxter integrable. 

In this integrable deformations framework ( |74| ) one notes that the 9x9 matrices H Jlyn+ i, considered for 
different values of their parameters V, J, t, commute all together and commute with P and P ■ H n ^ n+ i. 
Therefore one could imagine that the matrices R\ , and R2 , could belong to the same integrable algebraic 
variety. Is it possible to give a complete description of all the integrable deformations (|74|), and a 
complete description of the algebraic varieties associated with the Baxterisation of the largest integrable 
deformation of such quantum Hamiltonians ? These questions will be studied elsewhere. 



VIII. LET US BAXTERISE DIFFERENTIAL OPERATORS 



The Baxterisation is the building of the spectral parameters wherever they live : on an elliptic curve, 
on a rational surface, on a Jacobian variety, on an Abelian variety, ... The Baxterisation of R- matrices, 
or monodromy matrices, can always be performed systematically using various methods (analytical 
methods, by formal calculations, by visualization of the orbits pd| , ^2| , ...). Let us show, in this section, 
that these results can even be extended to infinite dimensional generalizations of monodromy matrices, 
namely when the block matrices in the i?-matrix are infinite dimensional linear operators (L-operators, 
local quantum Lax matrices ...). 



Let us consider, here, the L-operator corresponding to the Toda chain |5C 



L[X] 



X d 



(94) 



This can be seen as a generalization of the monodromy matrices ( ^T[ ) where the m x m matrices A, B, 
C and D become operators which can only have infinite dimensional representations. Let us introduce 
an "inverse" L-operator : 



I(L)[\>] 



(95) 



and let us perform the products of these two L-operators : 





L[\]-I(L)[\>] 
where : 



e- x 



N = 



A' + 4z 





"l N ' 




1 



(96) 



I(L)[X'].L[X] 
where : 



A' + £ j 



^ dx 


e x ' 




' 1 


" 


e~ x 







N' 


1 


d \ 
dx) 


(* 


d 

dx 


)■« 


—x 



(97) 



It is a straightforward calculation to see that the derivative operator ^ disappears in the off-diagonal 
N and N' operators, and that N and N' reduce to the operators of multiplication by the function 
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(A' + A — 1) • e x , and (A' + A — 1) • e x , respectively. The operators N and N' reduce to the null 
operator for A' = 1 — A. Therefore one has the following result : 



L[X] -I(L)[1 - A] = I(L)[1 - \].L[\] = 



1 
1 



(98) 



As it should the left "inverse" coincides with the right "inverse" . Let us use such inverse and combine 
it with the partial transposition t\, previously described. Let us now perform the partial transposition 
t\ which amounts to permuting the two off-diagonal operators : 



h(L)[X] 



e x 



(99) 



One easily finds, as a consequence of the nullity of two operators N\ and N2, when A' 

d 



-1 - A 



A-.,. = (\-±). e - x + e- x -(\' + ^- 
ax) \ ax 



and : N 2 = e x ■ A - 



dx 



+ (v 



d 

dx 



that the inverse of reads : 



I(h(L))[X'} = 



e~ x 

^ * + £ j 



where 



A' = -1 - A 



Let us perform, again, the partial transposition t\ on the previous L-operator : 

e x 



h(I(ti (L)))[X'] = 



(100) 



(101) 



The inverse of (101) reads 



I(ti(I(ti(L))))[X'' 



" ^ i 

e~ x 



(102) 



with A" = -A' +1. Actually : 

i 1 (/(ti(L)))[A']./(i 1 (I(t 1 (L))))[A" 



e x 
e~ x A' + 



dx 



\tf _d_ x 

A dx e 

e~ x 





1 0" 




_N 3 1_ 



e"* A' + -f 





" 1 iV 4 " 




1 



(103) 



where 



3 = e 



— X I \ll 



and : iV 4 = (V - -1 V e x + e x • ( X' + 

\ dx J \ dx 



d ) = (A" A' -I) ■<■• 



(104) 



Combining this result, and the previous result ( 100 ), one gets that I {t\{I (t\{L)))) has exactly the same 
form as L, but with A changed into A + 2 : 
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I(h(I(ti(L))))[X] = L[X + 2] 



(A +2)- £ e * 



(105) 



1 1\ 1 1\ is an infinite order transformation acting on the L-operator (p4|). This provides a first, and 
very simple, example of rational Baxterisation of the L-operator (p4|). The iteration of It\It\ yields a 
"trajectory" , in the space of the L-operator (94), corresponding to a straight line : A — > A + 2. It seems 
necessary, here, to perform (Iti) 2 , instead o f It %, to preserve the form of the L-operator. In fact one 
can only perform lt\ if one recognizes that ( |100| ) is like L(X + 1), up to a transformation c • s which 
commutes with the inversion / and with the partial transposition t\ : 



e" : 
e x -1-A 



dx _ 



(106) 



where c and s are the two following transformations (J7, V, W, T are any operator) 



and 



" u 


V ' 




' T 


w~ 


w 


T 


)- 


V 


u 



' u 


V ' 




' -u 


V 


w 


T 


)- 


w 


-T 



(107) 



Therefore It\ is associated with the shift A — ► A + 1. One sees that A actually plays the role of the 
spectral parameter for the Toda L-operator. 

Many other examples, corresponding to differential operators and L-operators, or local quantum Lax 
matrices associated with the XXX quantum Hamiltonian, yield similar calculations : for instance, a 



very simple example of L-operator, associated with the discrete self-trapping model (see for instance 51 1) 
can easily be Baxterised, the calculations being extremely similar to the one sketched in this section. 



A. Let us Baxterise higher derivatives 



Let us consider a (slightly) more complicated example of Baxterisation of differential operators. In the 
following we will denote by p the derivative operator -4- . More generally, let us consider the following 
L-operator : 



Ln — 



C{x) 



B(x) 




(108) 



where A^°\x) denotes the differential operator A (x) + Ai(x) -p H + A rn {x) -p m + ■ ■ ■ + An(x) -p N 

and where the two functions B(x) and C(x) are related by B(x) ■ C(x) = 1. Let us first perform the 
partial transposition t\ : 



h(L N ) 



A (x) + Ax(x)-p + 



+ A m (x)-p r ' 
B(x) 



+ ■•■ + A N {x)-p N C{x) 




(109) 



Let us try to find the inverse of this L-operator. Let us denote A^\x) another similar differential 
operator A (x) + Ai(x) ■ p + ■ • ■ + A m (x) ■ p m + ■ ■ ■ + Ajy(x) ■ p N . It is straightforward to see that 
t\(L]\r) ■ /(ti(Ljv)) is the identity matrix when /(ti(ijv)) reads: 



i(h(L N )) 





B(x) 



C(x) 



(110) 



with : 
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A (0 \x)-C(x) = C(x)-A (1) (x) 



(111) 



The fact that the left inverse and the right inverse identify is just a consequence of the fact that (111) 
is equivalent to B{x) ■ A^(x) — A^(x) • B(x) when B(x) ■ C(x) = 1. One easily sees that, up to 
the previous c and s transformations, the infinite order transformation I t\ reads : 



A^(x) B(x) 
C(x) 



<•( 



A^(x) B(x) 
C(x) 



(112) 



where A^(x) is deduced from A^(x) by (111). From now on let us note An, An, C the functions 
An{x), An{x) and C(x). We will denote C", C" and G^ N > the first, second and iV-th derivatives of 
C{x). For N = 3 the transformation A^(x) — » A^'(x) reads : 

C ■ A 3 = C ■ A 3 

C ■ A 2 = C ■ A 2 + 3 • C ■ A 3 

C ■ Ai = C ■ Ai + 2 • C ■ A 2 + 3 • C" ■ A 3 

C ■ A = C ■ A + C ■ At + C" ■ A 2 + C (3) • A 3 (113) 

Let us denote Si — C'/C, S2 = C"/C and S3 = C^'/C. One can straightforwardly associate to 
transformation ( |113| ) a 4 x 4 matrix : 

1 Si S2 S3 
1 2 Si 3S 2 
1 3 Si 
1 



M 



(114) 



The Baxterisation process corresponds to the iteration of this matrix. Note that this iteration yields a 
simple group structure : 

1 n- Si n ■ S2 + n ■ (n — 1) • Si 2 n ■ S3 + 3 • n ■ (n — 1) • Si S2 + n ■ (n — 1) • (n — 2) • Si 3 

1 2-n-Si 3-n- S 2 + 3-n- (n-l)- Si 2 
M n = (115) 

1 3-n-Si 

1 



If, inst ead of seeing 7(ti(ijv)) as an operator similar to Ljv, up to the previous c and s transfor- 
mations (107), one sees /(ii(Ljv)) as an operator similar to Ln, up to transformation c onl y, on e 
associates to It\ the 4x4 matrix —M, instead of M, yielding to a (—1)" factor in front of ( 115 ). 
If, similarly to the first Toda L-operato r exa mple, one just considers the iteration of (It\) 2 , one does 
not have a (—1)™ factor problem, but ( |115| ) is only valid for n even. Note that the "time reversal" 
I ■ t\ —tti-I corresponds to the same transformation as ( |113| ), but where the function C(x) is changed 

e r ' x the entries in the matrix 



into the function B(x) — 1/C(x). When C(x) is an exponential C(x) 
M are not functions of x but numbers. For N = 1 and C(x) = e x 
result associated to the 2x2 (shift) matrix : 



M 



1 1 
1 



one recovers the previous Toda 



(116) 



To sum up, transformation K 2 is associated with matri x M 2 (with M given by ( 114)). We are thus 
reduced to the analysis of the iteration of M given by (114). It is clear from (115) that the orbits of 
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K 2 in the "huge" functional space (Aq(x), A\[x), A2(x), ^3(2;)), correspond t o a r ational curve. This 



rational curve corresponds to the rational parametrization (Baxterisation) of (|115| ) which consists in 
considering the integer n as a real (or complex) number. 

Remark : It is worth noticing that all these Baxterisations of -L-operators have been performed 
without using any Yang-Baxter (RL L = L L R) hypothesis. These calculations are similar to the 
calculations on the sixteen vertex model where the iteration of K yields a foliation of the whole fifteen 
dimensional space of the model in elliptic curves, beyond the Yang-Baxter integrable cases |2q|. 



B. Some exercises of Baxterisation 



Of course all these calculations can be generalized to more general L-operators. We let the reader 
perform a few exercises of Baxterisation. 



• Let us Baxterise the L-operators : 



L = 



(117) 



where one recognizes ]52 ] the realization of the finite dimensional Lie algebra sl{2) in terms of differential 
operators : 



d_ 

dy 



h = -2y— + 21, f = -y 2 - — + 2ly 
ay ay 



(118) 



What is the "status" of the u and I parameters from a Baxterisation point of view : gauge parameters, 
spectral parameters ... ? 

• Show that the Baxterisation of the R-matrix 153] 

10 



R = 



u l-u 
10 
l-u u 



(119) 



yields & finite order group : the inverse I amounts to changing u into l/u, and transformation K 2 
t\ ■ I ■ t\ ■ I leaves the R- matrix (119) invariant : K 2 (R) = R. 



Show that the Baxterisation of the R-matrix 



R(q) 



yfq 

1 (g-l)/V« 

1 

^/q 



yields 



R(q, n) = 



jq 

1 {q-l)q n -^ 2 

1 

yfq 



(120) 



Hint : the inverse of R(q) is R(l/q) but K 2n (R(q)) = R(q,n) . 

• Recalling that the quantum group SL q (2) is a Hopf algebra with the antipode S given explicitly 
by [g : 





'A B~ 




D -B/q 


S(T) = S 






C D 




-qC A 



(121) 
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show that the antipodc S plays exactly the role of the inverse in the Baxterisation procedure. One will 
use the relations : 

CD-qDC=0, qB A — AB = 0, AD-DA= (q--)-BC, BC=CB (122) 

Hint: calculate the product T-S(T) and introduce the quantum determinant det g (T) = AD — qBC. 
This gives : 



T ■ S{T) ■- 
Introducing transformation t± : 



'A B~ 




D -B/q~ 




"1 






= det(T) • 

1 




C D 




-qC A 


1 



(123) 





'A 


B ' 




' A 


C~ 


ti{T) = h 












C 


D 




B 


D 



(124) 



let us consider the infinite dihedral group generated by these two involutions S and tx, or just the 
iteration of transformation t% ■ S. What is the "quantum group meaning" of this infinite dihedral group? 

• Let us recall the L-operator corresponding to the Liouville model on a lattice pa] : 



L = 



y/l -|- e 2Q+ih e P 



(125) 



where P and Q verify the Heisenberg commutation relation [P, Q] — ih, or the Hermann- Weyl 
relation e p — e lh e® e p . Try to Baxterise it, or remark that this problem reduces to the previous 
Hopf algebra, with a quantum determinant equal to 1 . 



• Let us recall the L-operator corresponding to the Sine-Gordon model (56| : 

u~ l/4m(«7A - Xv + 

l/4m (\v~ - v + /X) u + 



L 



(126) 



where it* and are the Weyl operators : — e ±l,9 / 4p and = e ±4 ' 5 / 2 ' ? . The spectral parameter 
A is explicit here. Show that the inverse transformation /, and the transformation K 2 , yield very 
simple multiplicative transformations on the spectral parameter A. In the classical limit this L-operator 
becomes: 



L = 



— 2 17 7T 

_ l/4m(Ae 1 /2#_ e -i/2i*/ x ) 



l/4m(Ae- 1 /2# _ e i/2 l 0/ A ) 
2 17 7T 



(127) 



In this classical limit what means the inversion relation / ? 

• Is it possible to Baxterise the universal R-matrix of U q {sl2) which corresponds to the formal se- 
ries i57t| : 



UU / — 1 \ fc 

R = Yl g fc(fc " 1)/2 ~L q 1 / 2H ^ H E k ^F k 



k=0 



[k}\ 



(128) 



31 



IX. CONCLUSION 



The Baxtcrisation procedure is very efficient, and powerful, and is not restricted to any Bethe Ansatz 
framework, or any particular "mathematical object" : one can actually Baxterise i?-matrices, mon- 
odromy matrices, L-operators, quantum Hamiltonians, ... in order to see the spectral parameter, and 
the integrability, become crystal clear. Typically, all the calculations sketched here, can be performed 
for R belonging to a quite general algebra R — X)i c » ' -^«> the generators Ai being not necessarily 
associated with semi-simple Lie algebras : they can, for instance, be elements of a Bose-Mesner alge- 
bra |38| associated with distance regular graphs ... The Baxterisation procedure corresponds to very 
simple calculations, namely performing the inverse of a matrix, or of some element of an algebra, and 
performing permutations of the entries of a matrix, and combining these two transformations to get 
a (generically) infinite order (birational, polynomial) transformation one studies as a discrete dynam- 
ical system. The Baxterisation procedure is a very powerful tool enabling to find, or simply analyse, 
Yang-Baxter integrable models. It is probably the quickest, and most powerful, way to get the solution 
(parametrization) of the Yang-Baxter equations, even if this parametrization is extremely complicated. 
Actually, the problem of the parametrization of the Yang-Baxter equations is often a quite difficult, or 
subtle, one. Recalling the free-fermion condition for the asymmetric eight vertex model, it is not always 
clear to see if a variable is a spectral parameter, an invariant (like the modulus of elliptic functions), 
or simply a gauge variable (5^^^] (see also A variable can be seen as a spectral parameter for a 

row-to-row transfer matrix, and as an invariant for the column-to-column transfer matrix, looking "like 
a gauge variable", but corresponding, in fact, to symmetries like ( |32| ) ... The Baxterisation procedure is 
a very efficient method to clarify these "subtleties" . 

It is striking to note that the Baxterisation procedure actually provides results beyond the Yang- 
Baxter integrable framework. The canonical elliptic parameterization of the sixteen vertex model is such 
a good example |2g| ] . What does the integrability of discrete symmetries of the parameter space of the 
model mean outside the Yang-Baxter integrable framework |62] ? Could it be possible that it could 
help to solve the model in the absence of Yang-Baxter integrability ? Is integrability restricted to Yang- 
Baxter integrability ? What is integrability? In our introduction we have recalled that the Yang-Baxter 
structure is a sufficient condition for the commutation of q N x q N transfer matrices (for any N) , which 
is a fundamental property for the integrability of the model. However, if one just wants to calculate the 
partition function of the model (largest eigenvalue), one does not need such commutation property in 
the whole q N dimensional space : if the transfer matrix can be block diagonalized, the commutation in 
some block, including the eigenvector corresponding to the largest eigenvalue, is sufficient to calculate 
exactly the partition function per site. The exact calculation of the partition function on the so-called 
"disorder solutions" f?3"| is such an example : the disorder solution "calculability" is not a Yang-Baxter 
integrability Jsl[ . Other examples exist in the literature called "quasi-integrability" : they correspond to 
Hamiltonians for which one can only find exactly the ground state and the corresponding eigenvalue. 

Along this line let us also recall the ideas developed by V. Jones on planar algebras |65| ] where he 
considers local relations, similar to the ABC = CBA Yang-Baxter relations, that could be sufficient 
to calculate global objects like generating functions equivalent to partition functions. He introduces, for 
instance, deformations of Yang-Baxter relations, like ABC = CBA + S, where S is "something" with 
"at most" two R matrices (see f55| ) . In such a framework one does not have a commutation of transfer 
matrices anymore. On could easily imagine to Baxterise the ABC = CBA + S relation in a similar 
way it can be done for the ABC = CBA Yang-Baxter relations ||[|. One could thus imagine the 
following situation : a parameterization of the deformed Yang-Baxter relation in terms of elliptic curves, 
no Yang-Baxter integrability stricto sensu, but, may be, a possibility to calculate exactly the partition 
function per site ! 

Leaving these speculative "dreams" let us just underline the fact that, beyond the Yang-Baxter frame- 
work, the Baxterisation procedure has already provided a very large number of new exact results for 
lattice models in statistical mechanics p6Lp7| an dj more generally, for discrete dynamical systems ptyXM . 
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X. APPENDIX A : POLYNOMIAL REPRESENTATIONS OF A — > M • A 



For N — 7, 11, one has a completely similar structure than the one depicted in section (VI), 
for the polynomial representation of the multiplication of the shift by N, namely (J x , J y , J z ) —* 

(4 N \ 4 N \ 4 N) ) where : 

j(N) _ j .p(N) (J j j \ 

J y ~ °V r y \ u xi Jy^ < J z) — 'Jy r x V J y> °zi "x) 

Jz ^ — Jz ' P% ^{Jxj Jyi Jz) — Jz ' J^x ^(Jz? Jx: Jy) (1^9) 



(7) 

For N — 7 polynomial Pi (J x , J y , J z ) reads 



■Pg* ^ ( J x i J y ; J z ) — 52 J y J 2 J j. 12 J y J g J r j. 495 J y J J ~\~ 260 J y J J ^ ~\~ 220 J y J ^ J ^ 

+188 Jl Jf Jf + 1680 J y 10 Jf Jl & - 42 Jf Jf Jl + 3192 jf jf J x 14 - 516 Jf Jf 4 

-5936 Jy 2 J z 18 J K 18 + 700 Jy 2 J 20 J K 16 - 364 j£ J 24 J K 18 - 66 J 4 J 20 J 24 + 12 J 2 J 22 J 24 - 25 20 Jy 4 J 22 J 12 

-2184 4° J 22 + 752 Jf Jf J 20 + 34 04 jf jf J 20 - 65 36 Jf J 14 J 20 + 5 7 04 J 14 jf jf - J 24 J 24 

+3864 J 14 Jf Jf + 1800 Jf J 24 Jf - 44 j£ J 22 J 2 f - J 24 J 24 + 7 J 24 J 24 + 756 J 8 J 22 Jf 

-1260 Jf J 24 Jf + 118 J 4 J 24 Jf + 168 Jf Jf Jf + 441 J 8 J 24 jf - 962 J 8 jf jf - 28 J 24 J 2 f J 2 

-1311 J 24 Jf Jl - 42 Jf Jf J 4 + 3756 J 20 jf jf 

+1800 J 24 Jf Jf + 484 J 24 Jf Jl - 1260 J 24 Jf Jf + 756 J 22 Jl Jf - 12 J 24 J 2 J 2 f - 28 J 22 J 24 J 2 X 
-2544 Jf Jf Jf - 9322 jf jf jf + 3756 jf jf jf + 12 jf J 2 jf 

-3928 Jf Jf Jf + 484 jf J 24 J% + 444 jf jf J s x + 792 jf jf J 24 + 744 J 22 jf jf - 2544 jf jf jf 

+444 Jf Jf Jl + 168 J 24 Jf Jf + 1734 jf jf J 8 - 52 J 22 jf J 4 + 48 Jf Jf Jf + 3864 jf jf jf 

+568 Jf Jf Jf + 3404 Jf Jf Jf + 744 jf J 22 jf - 5936 Jf Jf Jf + 752 jf jf jf 

+5704 Jj)8 Jf Jf - 516 Jf Jf Jl - 495 Jf Jl Jf + 220 Jf Jl Jf + 196 J 22 Jf J 4 + 824 jf jf jf 

+260 Jf Jl Jf - 1208 Jf Jf Jf - 924 jf jf jf + 792 jf jf jf - 1208 Jf Jf Jf 

+1680 Jf Jf Jf - 2184 Jf Jf Jf + 188 Jf J 6 Z Jf - 962 jf J s z jf + 824 jf jf jf 

-3928 Jf Jf Jf + 568 Jf Jf Jf + 700 Jf Jf Jf - 66 Jf Jf J 4 Z + 118 Jf Jf J 4 Z + 441 jf J 8 Z jf 

-2520 Jf Jf Jf - 1311 Jf Jf Jl + 3192 J 22 jf jf - 44 J 22 Jl jf - 364 J 24 Jl jf (130) 



The modulus ( [16[ ) is invariant by (129) corresponding to the multiplication of the shift by seven. 

The multiplication of the shift by eleven has a polynomial representation (J x , J y , J z ) — > 

(4 U \ Jy 11} , Jz} 11 ) satisfying Q for N = 11, where Pf 1] {J x , J y , J z ) is a homogeneous polyno- 
mial of degree 120 sum of 496 monomial expressions: 
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Pf 1] {J x , Jy, Jz) = -30045015 Jf Jf Jf + 2035800 Jf J 1 , 4 Jf + 142506 4° Jf Jf - 27405 Jy Jf Jf 
+54627300 jf jf jf - 435 jf J 4 jf + 54627300 jf jf jf + ■ ■ • (131) 

The actual expression of Px\j x , Jy, Jz) will be given elsewhere. 

When N is not a prime number one has slightly modified results: one does not have ( |129| ) anymore. 
Actually, the polynomial representation of the multiplication of the shift by six can be obtained in many 
different ways, namely substituting in j( 3 \ or in jM, or by various eliminations between 
various biquadratic curves I\. The result reads {J x , J y , J z ) — > {jf\ Jy & \ Jz 6 " 1 ) where : 

Jx = Jx ' \ Jx, Jy, Jz) 

Jy ^ = Jy ^ ^ {Jx, Jy, Jz) Jy ^ P x {Jy, Jz, Jx) 

Jz ^ = Jz ^ J'z ^{Jx, Jy, Jz) — Jz^ ' Px ^{Jz, Jx, Jy) (132) 
SUCh that jf\j x , Jy, Jz) = Jx e) (Jy, Jz, Jx) and jf\j x , Jy, J z ) = ji & \j z , Jx, Jy)- The P^'s 

and the J^'s are functions of J x , Jy and Jf One also has J X {J X , J y , J z ) = Jx & \j x , Jz, Jy), 

Jy {Jx, Jy, Jz) = Jy {Jz, Jy, Jx), Jz {Jx, Jy, Jz) — Jz {Jy, Jx, Jz) and P x {Jx, Jy, Jz) — 

P^IT J J \ P (6 V 7 7 T ) — P^(T J J] P^tJ J J) — P^tJ J J) 
r x \ ,J x, Jz, Jy) , r y \Jx, Jy, Jz) — r y \Jz, Jy, Jx) , r z \Jx, Jy, Jz) — r z \Jy, Jx, Jz)- 

All these expressions can thus be deduced from P x {J x , Jy, Jz) ■ 
J'x ^ {Jx, Jy, Jz) — 70 J z J x Jy — 250 J z J x Jy — 8 J z J x Jy + 8 J z J x Jy +184 J z J x J y 

I 2 8 J % x y ^ ^ z J x y ^ ^ z J x y ^ z x ,J y ^ z x y ^ z x y ^ z x 

+228 J z J x Jy -+- 40 J z J x Jy 4" 70 Jy J z J% — 56 Jy J z J x + 88 Jy® J\ J x + 28 Jy ^ e/^ t/^. + 56 <^y^ ^ J x 
-40 J* 4 J x 4 + 56 Jf Jf J% + 88 Jf Jf Jl - 8 Jf Jf Jl + 28 Jf Jf J 4 - 56 Jf Jf J* x 
+40 jf jf J z + 184 jf jf J z + 228 jf jf J z — 72 jf jf J® — 68 Jf Jf J z — 88 jf Jf J z 

56 jy J x j z + 40 jy j x j z + 8 jy (^.J, 56 jy j^. £ + 28 jy 88 j z j rj. jy 316 jy j^. 

— 8 jy j x j z + 144 j z j x jy — 400 j z J x jy ^ l ^ j g j^. j y ~\~ J ^ ^ + j y j^. ~ 8 ^ j g ^ 

Note that J^ 6 " 1 is a homogeneous polynomial expression of degree 36. 

Similarly, the polynomial representation of the multiplication of the shift by nine rGcids {Jx-> Jyi Jz) ^ 
(Ji 9) ,4 9) ,ji 9) ) where: 

J x ^ = J x P x {J X , Jy, Jz) 

Jy ^ = Jy ^ Py (^li Jy, Jz) — ^x i^V Jz, Jx) 

J z ^ = J z P z {J Xi Jy, J z ) — Jz Px {J%1 Jx, Jy) (133) 

The expression of Px 9 \jx, Jy, Jz) will be given elsewhere. It is, again, a function of Jj, Jy and J z . It 

is the sum of 190 monomial expressions of degree 72, jf^ being a homogeneous polynomial expression 
of degree 81. 

XI. APPENDIX B : FINITE ORDER CONDITIONS 

It can be shown that the points of the Baxter model on the algebraic varieties : 
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^ ^ ( ^sc ; Jyi J z) — J x Jy H~ J z J x ~^ J z Jy •> Vz ( ) ^Vi ^ z) — Jx Jy H~ J z J x J z Jy 

are actually such that K G (R) = ( • R. One has the following factorization property : 

V^(4 2 \ jf , jf ) = (134) 

V&{J X , Jy, J z ) ■ V®\J X , Jy, Jz) ■ V^(J X , Jy, Jz) ■ V^(J X , Jy, Jz) = (135) 

where 14 ; (J X , J a , J 2 ) = V z '(J y , J z , J x ) and Vjj '(J x , J y , J z ) = Vz (J z , Jx, Jy)- One has the 
relation : 

^(jf , jf , jf ) - Pf(J„ Jy, J z ) = (136) 



Note that the points of the Baxter model on the algebraic varieties P x (J x , Jy, Jz) = 0, or 
Py 3) (J x , J v , Jz) = 0, or pJ 3) (J x , J y , J z ) = 0, are such that K 12 {R) = ( ■ R. Relation jl36| ) is in 

(2) (2) (2} 

agreement with the fact that ( J x , J y , J z ) — * (Jx , Jy , Jz ) is a repr esentation of the shift doubling. 
The points of order six (namely K 6 (R) = ( ■ R) correspond to (|l 34j) , their image by the shift dou- 



bling (J x , J y , J z ) -> (J x 2 >/f , jf ) giving P^(J X , J y , J z )-Pf(J,, J y , J Z )-P^\J X , J y , Jz) = 0, 



together, of course, with (134) 



Let us note that the points of the algebraic variety V z (J x , J y , J z ) — are (projectively) of order 
three for the Baxter model : K 3 (R) = £ • R. The points of the Baxter model on the variety : 

V^( J T T } — T 2 J T 3 — 9 T 2 T 2 T 2 4- T 2 T 3 J 4- J J 2 T 3 — J J 3 T 2 — J 3 T 3 — T 3 T J 2 

'y \»X, J y, <Jz) — J z J x J y 4 o z J x J y ' ° z J x J V ' " z J x J y " z J x J y " z J y ° z J x J y 

~^~Jz Jx Jy + J z Jx ~ Jx Jy = (137) 



are of order five K a (R) — ( ■ R. The points of the Baxter model on the algebraic variety : 

Vy ^ {Jx, Jy, Jz) — Vy ^ {Jy, Jx, Jz) = J z Jx Jy ~ 2 J z J x Jy + J z J x Jy + J z J x Jy 

— J Z J x Jy — J z J X ~ J Z Jx Jy J Z Jx Jy + J ' z ^ J/ _ ^ X J J/ = (1^8) 

are of order ten : K 10 {R) = ( ■ R. The point of the Baxter model on the variety Px (Jx, Jy, Jz) = 0, 
Pip (Jx, Jy, Jz) = 0, or Pz 5 \j x , J y , J z ) = 0, are of order twenty : K 20 (R) = £ ■ R. Let us note that: 

vf> (jf , jf , jf ) + p^(j x , j y , j z ) = o, vf> (jf , jf , Jf ) + pf (j XJ j„ j,) = o (139) 

(2) (2) (2} 

which is in agreement with the fact that (J x , J y , J z ) — > (Jr , Jy , Jz ) represents the shift doubling. 



Remark : Relation (139) is in agreement with the shift doubling, however one seems to have an 



apparent contradiction with previous relations. The points of (137), namely Vy (J x , J y , Jz) = 



being of order five, one expects that their image by the shift doubling will give points of order ten, 



like (138), and not points of order twenty, like Py 5 \j x , J y , J z ) — 0. In fact, this algebraic variety is 



an order five algebraic variety only when restricted to the Baxter model. For the sixteen vertex model 
one can actually show that Vy 5 \j x , J y , J z ) = and Vy 5 \j x , J y , J z ) = 0, are (codimension-one) 
algebraic varieties of order ten, the algebraic varieties of order five being higher codimension algebraic 
varieties. 
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